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Managementuittreksel 


Titei  Berekening  van  het  in-  en  extern  risico  van  explosievenopslag  met 

behulp  van  ‘RISKANAL’ 

Auteur(s)  Ir.  P.G.J.  Timmers 

Datum  :  juli  1998 

Opdrachtnr.  A97KL474 

Rapportnr.  PML  1998-A23 


Het  computerprogramma  ‘RISKANAL’  is  door  het  TNO  Prins  Maurits  Laboratori- 
um  in  opdracht  van  het  Ministerie  van  Defensie  ontwikkeld  om  risico ’s  verbonden 
aan  de  opslag  van  munitie  en  explosieven  te  kwantificeren.  Via  een  interactieve 
vragenlijst  wordt  de  lay-out  van  een  munitiemagazijnencomplex  in  het  programma 
ingevoerd  en  kan  het  interne  en  exteme  risico  verbonden  aan  de  opslag  berekend 
worden.  Voor  een  nieuwe  opslag  is  het  tevens  mogelijk  om  een  optimale  lay-out  te 
bepalen,  wat  resulteert  in  een  zo  laag  mogelijk  risico.  Ook  kan  het  effect  van  een 
andere  belegging  van  de  explosieven  op  het  extern  risico  direct  in  kaart  worden 
gebracht. 

Op  verzoek  van  het  Ministerie  van  Defensie,  DMKL/sectie  Milan,  bureau  AMV 
werd  een  nederlandstalige  beschrijving  gegeven  van  de  werking  van  het  program¬ 
ma  ‘RISKANAL’.  Hiermee  wordt  een  duidelijk  overzicht  gegeven  van  de  wijze 
van  aanpak  en  de  methodes  die  door  het  progranuna  gebruikt  worden.  Het  rapport 
kan  tevens  gebruikt  worden  als  handleiding  voor  het  gebruik  van  ‘RISKANAL’  of 
voor  de  interpretatie  van  de  resultaten. 

Het  programma  berekent  het  intern  en  extern  risico  van  een  explosievenopslag. 

Met  intern  risico  wordt  de  mogelijkheid  van  een  initiatie  in  een  magazijn  als  ge- 
volg  van  een  ontploffend,  nabij  gelegen  magazijn  bedoeld.  Het  extern  risico  geeft 
het  individueel  en  groepsrisico  voor  personen  in  de  omgeving  als  gevolg  van  een 
exploderend  magazijn.  Het  programma  is  gebaseerd  op  een  combinatie  van  de 
aanbevelingen  van  de  NATO-werkgroep  AC/258  en  de  letselcriteria  zoals  beschre- 
ven  in  het  Nederlandse  ‘Groene  Boek’.  De  AC/258-aanbevelingen  zijn  gebruikt 
voor  de  berekening  van  de  explosie-effecten  en  voor  de  berekening  van  het  intern 
risico.  Voor  de  berekening  van  de  fragment-  en  brokstukdichtheid  wordt  gebruik- 
gemaakt  van  door  het  TNO  Prins  Maurits  Laboratorium  afgeleide  formules,  geba¬ 
seerd  op  grootschalige  testen  in  USA  (ESKIMO,  Yuma).  De  formules  uit  het 
groene  boek  worden  gebruikt  voor  de  berekening  van  de  gevolgen  voor  personen  in 
de  omgeving.  Vervolgens  kan  met  behulp  van  de  kans  op  een  ongeval  het  extern 
risico  bepaald  worden. 


Uit  de  beschrijving  van  het  programma  worden  de  volgende  conclusies  getrokken. 

•  Door  het  interactieve  karakter  van  het  programma  kunnen  de  input-gegevens 
eenvoudig  worden  veranderd  en  kunnen  de  effecten  van  deze  veranderingen 
snel  zichtbaar  gemaakt  worden.  Hierdoor  is  het  mogelijk  om  het  effect  als  ge- 
volg  van  wijzigingen  in  de  lay-out  van  het  munitiemagazijnencomplex  in  kaart 
te  brengen  en  zo  te  komen  tot  een  optimaal  ontwerp  met  betrekking  tot  veilig- 
heid. 

•  Het  programma  is  zodanig  opgezet  dat  de  gebruikte  formules/relaties  door 
ingewijden  in  het  programma  vrij  eenvoudig  kunnen  worden  vervangen  door 
nieuwe  relaties  verkregen  uit  nieuwe  inzichten.  Men  kan  hierbij  denken  aan 
nieuwe  inzichten  met  betrekking  tot  schadecriteria,  brokstukkenuitworp,  initia- 
tiekans,  enzovoort. 

•  De  absolute  getallen  in  de  output-file  dienen  met  enige  voorzichtigheid  te 
worden  geinterpreteerd.  Doordat  wereldwijd  slechts  beperkte  experimentele  ge- 
gevens  beschikbaar  zijn,  zijn  veel  van  de  gebruikte  formules  slechts  beperkt  ge- 
valideerd.  De  onnauwkeurigheid  van  de  resultaten  kan  door  de  gebruikte  for¬ 
mules  relatief  groot  zijn.  Om  dit  te  verbeteren  is  experimenteel  onderzoek 
noodzakelijk.  De  verkregen  kwantitatieve  resultaten  zijn  echter  een  goed  hulp- 
middel  om  waarden  voor  het  extern  risico  relatief  met  eUcaar  te  vergelijken. 
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Samenvatting 


Sinds  het  begin  van  de  jaren  ’80  wordt  voor  risicoanalyse  van  explosievenopslag 
door  het  TNO  Prins  Maurits  Laboratorium  (TNO-PML)  gebruikgemaakt  van  het 
computerprogramma  ‘RISKANAL’.  Dit  programma  werd  door  TNO-PML  zelf 
ontwikkeld  en  steeds  aangepast  aan  de  meest  recente  kennis  met  betrekking  tot  de 
risicoanalyse.  In  dit  rapport  wordt  een  beschrijving  gegeven  van  het  programma 
‘RISKANAL’,  de  opbouw  en  gebruikte  rekenmethodes,  waardoor  het  rapport 
tevens  gebruikt  kan  worden  als  handleiding  voor  het  gebruik  van  het  programma. 
‘RISKANAL’  berekent  de  interne  en  exteme  risico’s  behorende  bij  een  explosie- 
venoplag.  Met  intern  risico  wordt  bedoeld  de  mogelijkheid  dat  een  magazijn  wordt 
geinitieerd  door  een  ontploffend  emaast  gelegen  magazijn.  Het  extern  risico  geeft 
het  aantal  slachtoffers  in  de  omgeving  als  gevolg  van  een  exploderend  magazijn. 
Het  programma  gebruikt  de  veiligheidsafstanden  volgens  de  aanbevelingen  van  de 
NATO-werkgroep  AC/258  [2]  voor  de  berekening  van  het  intern  risico  en  de 
probit-functies  zoals  gegeven  in  het  ‘Groene  Boek’  [3]  voor  de  berekening  van  de 
gevolgen. 

Het  uiteindelijke  resultaat  van  het  programma  is  het  gemiddelde  aantal  slachtoffers 
(letaliteiten)  dat  verwacht  kan  worden  wanneer  een  magazijn  op  een  complex  zou 
exploderen.  Met  deze  getallen  kan  het  groepsrisico  worden  berekend  door  de  kans 
op  een  accidentele  explosie  erbij  te  betrekken.  Op  grond  van  deze  informatie  kan 
beslist  worden  of  een  situatie  al  of  niet  aanvaardbaar  is. 
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1  Inleiding 


Sinds  het  begin  van  de  jaren  ’80  wordt  voor  risicoanalyse  van  explosieven-/ 
munitieopslag  door  het  TNO  Prins  Maurits  Laboratorium  (TNO-PML)  gebruikge- 
maakt  van  het  computerprogramma  ‘RISKANAL’.  Het  programma  ‘RISKANAL’ 
werd  door  TNO-PML  ontwikkeld  om  risico’s  verbonden  aan  de  opslag  van  munitie 
en  explosieven  te  kwantificeren.  Hiermee  kunnen  interne  en  exteme  risico’s  bij 
bestaande  opslag  van  explosieven  en  munitie  berekend  worden  en  is  het  tevens 
mogelijk  voor  een  nieuwe  opslaglocatie  een  optimaal  ontwerp  te  realiseren,  wat 
resulteert  in  een  zo  laag  mogelijk  risico.  Daamaast  kan  het  effect  van  een  andere 
belegging  op  het  extern  risico  direct  in  kaart  worden  gebracht,  waardoor  de  moge- 
lijkheid  onderzocht  kan  worden  om  door  aanpassing  van  de  belegging  het  extern 
risico  te  verlagen. 

Door  het  Ministerie  van  Defensie,  DMKL/sectie  Milan,  bureau  AMV  werd  ge- 
vraagd  een  nederlandstalige  beschrijving  te  geven  van  de  werking  van  het  pro- 
gramma  ‘RISKANAL’.  Hiermee  wordt  een  duidelijk  overzicht  gegeven  van  de 
wijze  van  aanpak  en  de  methodes  die  door  het  programma  gebruikt  worden,  Het 
rapport  kan  tevens  gebruikt  worden  als  handleiding  voor  het  gebruik  van 
‘RISKANAL’  of  voor  de  interpretatie  van  de  resultaten. 

Het  rapport  zal  eerst  een  algemene  beschrijving  geven  van  de  werking  van  het 
programma  RISKANAL.  Hiema  zullen  de  verschillende  onderdelen  van  het  pro¬ 
gramma  elk  afzonderlijk  behandeld  worden. 
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2  Het  programma  ‘RISKANAL’ 


Het  programma  ‘RISKANAL’  is  ontwikkeld  om  de  interne  en  exteme  risico’s  te 
berekenen  bij  munitie-  of  explosievenopslag.  Met  intern  risico  wordt  de  mogelijk- 
heid  tot  sympatische  detonaties  bedoeld  indien  een  magazijn  op  een  magazijnen- 
complex  explodeert.  Bij  een  sympatische  detonatie  vindt  initiatie  van  de  in  een 
magazijn  opgeslagen  explosieven  plaats  door  het  exploderende  naastgelegen  maga¬ 
zijn.  Het  exteme  risico  is  gedefinieerd  als  de  kans  op  letale  verwondingen  voor 
personen  in  de  directe  omgeving  van  de  opslaglocatie  per  jaar.  Er  wordt  een  indivi- 
dueel  risico  en  een  groepsrisico  berekend. 

Voor  de  berekening  van  het  interne  risico  en  het  kwantificeren  van  de  explosie- 
effecten  (‘blast’  en  hittestraling)  wordt  gebmikgemaakt  van  de  ACy258- 
aanbevelingen  [2].  Het  kwantificeren  van  de  uitworp  van  fragmenten  en  brokstuk- 
ken  wordt  gedaan  met  behulp  van  door  TNO-PML  afgeleide  empirische  formules. 
Voor  de  bepaling  van  de  gevolgen  (extern  risico)  worden  de  probit-fimcties  voor 
letaal  en  emstig  letsel  uit  het  Nederlandse  ‘Groene  Boek’  [3]  gebmikt.  Een  probit- 
functie  geeft  de  kans  op  letsel  als  functie  van  de,  voor  dat  letsel,  meest  belangrijke 
parameters.  Zo  zijn  de  meest  belangrijke  parameters  voor  letale  longschade  de 
bereikte  piekoverdmk  en  impuls  en  kan  uit  de  probit-functie  voor  letale  longschade 
het  percentage  sterfgevallen  worden  afgeleid  bij  een  bepaalde  combinatie  van 
piekoverdmk  en  impuls. 

De  belangrijkste  aannamen  voor  de  berekening  van  de  explosie-effecten  zijn  als 
volgt. 

•  ‘Blast’-parameters  zijn  afgeleid  van  halfbolvormige  TNT-ladingen  die  op 
maaiveldniveau  en  in  de  open  lucht  tot  ontploffing  worden  gebracht. 

•  Effectberekeningen  gaan  uit  van  een  vlak  terrein.  De  invloed  van  hellingen, 
vegetatie  en  andere  obstakels  wordt  niet  meegenomen. 

•  De  ‘blast’-parameters  zijn  geldig  voor  geschaalde  afstanden  van  0,0674  tot 
40,0  m/kg  De  geschaalde  afstand  legt  een  relatie  tussen  de  hoeveelheid  ex- 
plosieve  stof  en  de  afstand  waarbij  een  bepaalde  piekoverdmk  optreedt.^  Deze 
relatie  is  weergegeven  in  tabel  1. 

Met  de  formules  uit  ACy258,  part  II  [2]  wordt  door  het  programma  voor  elk  do- 
normagazijn  met  zijn  totale  betrokken  netto  hoeveelheid  explosieven  de  invallende 
en  gereflecteerde  dmk,  de  invallende  en  gereflecteerde  impuls,  de  duur  van  de 
positieve  fase,  de  stralingsintensiteit  en  de  stralingsduur  berekend.  De  fragment- 
dichtheid  wordt  berekend  met  empirische  formules  afgeleid  van  experimentele  data 
verkregen  tijdens  de  ESKIMO  [4,  5]  en  Yuma  [6]  tests  (10.900  kg  M107, 


1  Uit  de  natuurwetten  blijkt  dat  als  de  afstand  tot  de  detonatie  wordt  verdubbeld  er 
8  keer  zoveel  explosief  nodig  is  ora  daar  dezelfde  schokgolfgrootte  te  krijgen.  De 
geschaalde  afstand  is  daarom  gedefinieerd  als  SD  =  m/q*^^  (met  ra  is  de  afstand  in 
meter  en  q  is  de  hoeveelheid  explosieve  stof  in  kg).  Met  andere  woorden  een  ex- 
plosie  van  1000  kg  explosief  op  100  meter  afstand  geeft  dezelfde  piekoverdruk  als 
een  explosie  van  10  kg  op  21,5  meter  afstand. 
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10.900  kg  Ml  17  en  2970  kg  Ml  17).  Met  de  verkregen  resultaten  vindt  een  schat- 
ting  plaats  van  de  effecten  op  omliggende  objecten,  waarbij  gebruikgemaakt  wordt 
van  de  probit-functies  uit  bet  Nederlandse  ‘Groene  Boek’  [3]. 

Tabel  1:  Relatie  tussen  geschaalde  af stand  (SD)en  de  te  verwachten  piekoverdruk 

volgens  NATO  ACI258-aanbevelingen  [2]. 


Geschaalde  afstand  (SD) 

(Q  in  kg,  afstand  in  m) 
(m/kg^^3) 

Venwachte  piekoverdruk  (side-on) 

(bar) 

55,5 

0,015 

44,4  tot  33,3 

0,02  tot  0,03 

22,2 

0,05 

14,8 

0,09 

9,6 

0,16 

8,0 

0,21 

7,2 

0,24 

3,6 

0,70 

2,4 

1,80 

Het  programma  ‘RISKANAL’  is  opgedeeld  in  vier  delen.  In  het  eerste  deel  van  bet 
programma  dient  de  situatie  gescbetst  te  worden.  Dit  betekent  dat  potentieel  ge- 
vaarlijke  opslagfaciliteiten  (donors)  en  omliggende  objecten/personen  welke  mo- 
gelijk  scbade  of  letsel  kunnen  ondervinden  van  de  exploderende  donor  (acceptors) 
gedefinieerd  moeten  worden.  In  bet  tweede  deel  van  bet  programma  wordt  bet 
‘intern  risico’  berekend.  Dit  betekent  dat  wordt  nagegaan  of,  en  zo  ja  welke,  donors 
een  sympatiscbe  detonatie  veroorzaken.  De  boeveelbeid  bij  een  detonatie  betrok- 
ken  explosief  wordt  aan  de  band  van  deze  resultaten  aangepast.  In  bet  volgende, 
derde,  deel  van  bet  programma  wordt  bet  ‘extern  risico’  berekend.  Dit  betekent  bet 
kwantificeren  van  de  explosie-effecten  van  elke  donor,  de  gevolgen  voor  elke 
acceptor  en  bet  bepalen  van  een  gemiddelde  kans  op  letaliteit.  Tenslotte  wordt  bet 
gebeel  weergegeven  in  een  output-file  waarin  de  risico ’s  worden  aangegeven. 


Om  tot  een  goede  bescbrijving  van  bet  programma  te  komen  is  bet  waarscbijnlijk 
bet  eenvoudigst  om  de  lijn  van  het  programma  te  volgen.  In  het  vervolg  van  dit 
hoofdstuk  zullen  de  vier  onderdelen  van  het  programma  daarom  afzonderlijk  aan 
bod  komen  en  uitvoerig  worden  behandeld. 


2.1  Definitie  van  donors  en  acceptors 

In  dit  gedeelte  van  het  programma  wordt  door  middel  van  het  stellen  van  vragen 
een  beeld  gescbetst  van  de  situatie  op  en  om  het  explosieven-  of  munitieopslagter- 
rein.  Alle  stappen  zullen  nu  worden  besproken. 

Om  snel  een  indruk  te  krijgen  van  de  globale  effecten,  zonder  eerst  alle  gegevens 
in  te  hoeven  voeren,  start  het  programma  met  de  mogelijkheid  om  voor  maximaal 
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vijf  verschillende  netto  hoeveelheden  explosieven  (gevarenklasse  1.1)  het  effect 
(kans  op  letaliteit/gewonden)  uit  te  zetten  als  functie  van  de  afstand.  Indien  deze 
grafieken  gewenst  zijn  dan  zal  het  programma  door  middel  van  een  subroutine  de 
situatie  ‘een  magazijn  en  e6n  persoon  in  een  huis  of  in  het  vrije  veld’  berekenen 
voor  verschillende  afstanden.  Dit  kan  worden  uitgevoerd  met  of  zonder  het  effect 
van  fragmenten  en  brokstukken.  In  bijlage  A  worden  enkele  voorbeelden  van  deze 
figuren  gegeven.  Indien  verder  geen  figuren  meer  gewenst  zijn,  kan  het  programma 
hiema  worden  afgesloten. 

Als  het  programma  niet  wordt  afgesloten  wordt  vervolgens  de  naam  van  de  input- 
file  gevraagd.  Deze  heeft  maximaal  negen  karakters  en  een  extensie  .DAT.  Als  de 
input-file  al  bestaat  wordt  deze  gelezen.  Bij  een  nieuwe  input-file  zal  deze  worden 
gecreeerd. 

Het  zal  duidelijk  zijn  dat  een  grote  hoeveelheid  data  nodig  is  om  een  munitiemaga- 
zijnencomplex  (MMC)  en  zijn  omgeving  zo  te  kunnen  beschrijven  dat  de  output 
realistisch  is.  Aangezien  deze  data  voor  de  meeste  opslagterreinen  in  Nederland  al 
in  het  juiste  format  beschikbaar  zijn  is  geen  bijzondere  aandacht  besteed  aan  een 
programmaonderdeel  dat,  bijvoorbeeld  met  een  grafisch  tablet,  het  maken  van  zo’n 
input-file  vergemakkelijkt. 

De  input-file  dient  informatie  te  bevatten  over  het  complex  zelf  (type,  orientatie  en 
positie  van  de  magazijnen  met  zijn  inhoud  aan  explosieve  stoffen  en  de  geva- 
ren(sub)klasse  daarvan)  en  over  de  omgeving  (type  en  positie  van  de  objecten  en 
het  aantal  personen  per  object).  De  afmetingen  en  de  orientatie  van  een  magazijn 
(donor)  worden  gedefinieerd  door  zijn  coordinaten.  Een  magazijn  wordt  bepaald 
door  drie  codrdinaten.  De  eerste  twee  coordinaten  bepalen  de  wand  met  de  ingang 
van  het  magazijn,  de  derde  coordinaat  is  voor  de  achterwand.  Deze  coordinaten 
zijn  nodig  om  de  orientatie  vast  te  leggen,  uitgaande  van  een  rechthoekig  magazijn. 
Deze  ligging  is  zeer  belangrijk,  aangezien  de  wand  met  de  ingang  het  meest  kwets- 
baar  is,  zowel  voor  exteme  effecten  als  in  het  geval  van  een  interne  explosie. 

Elk  object  (acceptor)  in  de  omgeving,  behalve  een  weg,  wordt  vastgelegd  met  een 
paar  coordinaten.  Wegen  worden  gedefinieerd  door  het  begin-  en  eindpunt  van 
rechte  lijnen  (minimaal  25  meter).  Een  weg  met  drie  rechte  stukken  kan  worden 
beschreven  met  vier  coordinaten.  Voor  een  kromming  in  een  weg  dient  een  ver- 
eenvoudiging  plaats  te  vinden  door  deze  op  te  delen  in  rechte  stukken  van  maxi¬ 
maal  25  meter  lengte. 

Als  een  nieuwe  input-file  gemaakt  moet  worden,  vraagt  het  programma  interactief 
om  de  volgende  gegevens: 

•  het  aantal  magazijnen; 

•  het  aantal  omringende  objecten; 

•  het  totale  aantal  objectcoordinaten; 

•  het  gemiddelde  aantal  personen  per  unit; 

•  de  naam  van  het  complex; 

•  voor  elk  magazijn: 
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-  het  type  van  het  magazijn; 

-  de  naam  van  het  magazijn; 

-  de  explosieve  inhoud  van  het  magazijn; 

-  de  gevarenklasse  van  de  inhoud; 

-  drie  hoekcoordinaten  van  het  rechthoekige  magazijn; 

•  voor  elk  omringend  object: 

-  de  naam  van  het  object; 

-  het  type  object; 

-  het  nummer  in  ‘array’  van  de  eerste  coordinaten  (k)^; 

-  het  nummer  in  ‘array’  van  de  laatste  coordinaten  (n)^; 

-  het  aantal  units  voor  het  object  (aantal  personen)  of  de  24-uurs  verkeersin- 
tensiteit  van  auto’s  en  fietsers/voetgangers  voor  wegen; 

-  de  centrumcodrdinaten  van  het  object  of  de  eerste  coordinaten  van  een  weg; 

-  de  volgende  (n-k)  coordinaten  van  een  weg. 

Voor  het  type  magazijn  is  een  aantal  mogelijkheden  in  het  programma  opgenomen, 
elk  met  een  karakteristieke  code: 

17  :  igloos  ontworpen  voor  7  bar; 

13  :  igloos  ontworpen  voor  3  bar; 

ID  :  met  aarde  overdekte  constructie  met  een  frontpui  en  deuren  die,  indien  ze 
gericht  zijn  naar  een  potentieel  explosiegevaarlijke  locatie,  bestand  zijn  te- 
gen  hogesnelheidsprojectielen; 

IB  :  met  aarde  overdekte  constructie  met  een  deur  die,  indien  deze  gericht  is  naar 
een  potentieel  explosiegevaarlijke  locatie,  gebarricadeerd  is; 

IE  :  niet  gebarricadeerde  met  aarde  overdekte  constructie; 

BR  :  zware  constructies  (0,7  m  beton,  steen  of  equivalent)  met  een  beschermend 
dak  (0,15  m  beton  met  bijbehorende  draagconstructie)  en  een  gebarricadeer¬ 
de  deur  indien  deze  gericht  zijn  naar  een  potentieel  explosiegevaarlijke  loca¬ 
tie; 

BD  :  zware  constructie  zonder  beschermend  dak  maar  met  een  gebarricadeerde 
deur  indien  gericht  naar  een  potentieel  explosiegevaarlijke  locatie; 

OB  :  aan  drie  zijden  gebarricadeerde  openluchtopslag  of  lichte  constructie; 

OT  :  aan  vier  zijden  gebarricadeerde  openluchtopslag  of  lichte  constructie; 

OP  :  niet  gebarricadeerde  openluchtopslag  of  lichte  constructie. 

Ook  voor  het  invoeren  van  omringende  objecten  is  een  aantal  typen  mogelijk,  dat 
elk  met  een  karakteristieke  code  wordt  aangeduid: 

HS  :  bewoond  huis; 

HF  :  bewoond  huis  met  meer  dan  vier  verdiepingen; 

HU  :  personen  in  tenten,  parken,  enzovoort  (onbeschermde  (‘unprotected’)  perso¬ 
nen); 

RD  :  openbare  verkeersweg; 

PL  :  POL-faciliteiten  (Petroleum,  Oils,  Lubricants). 


2  Vcwr  alle  objecten,  behalve  wegen,  geldt  dat  n  gelijk  is  aan  k. 
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Met  een  vooraf  bepaald  gemiddeld  aantal  personen  per  unit  kan  het  totale  aantal 
personen  per  object  worden  bepaald  door  een  aantal  units  te  koppelen  aan  elk 
omringend  object.  Als  bijvoorbeeld  het  gemiddeld  aantal  personen  per  unit  2,5 
bedraagt,  kan  een  groep  van  10  personen  in  het  park  worden  gedefmieerd  door 
4  units  (4  HU).  Dit  betekent  wel  dat  gegevens  beschikbaar  moeten  zijn  over  het 
gemiddeld  aantal  personen  per  huis  of  in  een  park.  Voor  een  openbare  weg  dient  de 
24-uurs  verkeersintensiteit  voor  auto’s  en  voor  fietsers/voetgangers  bekend  te  zijn 
of  te  worden  geschat.  Het  programma  neemt  voor  auto’s  en  vrachtverkeer  een 
gemiddelde  snelheid  aan  van  60  km/uur,  en  voor  fietsers  15  km/uur. 

Bij  het  beschrijven  van  de  omgeving  is  het  van  belang  te  weten  tot  welke  afstand 
van  de  donors  de  omgeving  in  een  berekening  meegenomen  moet  worden.  Deze 
afstand  is  namelijk  afhankelijk  van  de  explosie-effecten  en  dus  van  de  geva- 
ren(sub)klasse  en  hoeveelheid  van  het  opgeslagen  explosief.  Voor  een  opgeslagen 
hoeveelheid  1.1  explosief  geldt  een  zonering  volgens  de  zogenaamde  nota  van  van 
Houwelingen  [7].  Hierin  wordt  een  indeling  gemaakt  naar  A-zone,  B-zone  en 
C-zone,  afhankelijk  van  de  afstand  tot  de  donor.  Buiten  de  grootste  zone,  de 
C-zone,  zijn  geen  restricties  meer  geldig.  Het  lijkt  daarom  aannemelijk  om  deze 
zone  als  grens  te  nemen  voor  de  input-file  van  RISKANAL.  De  C-zone  ligt  op 
tweemaal  de  afstand  voor  huizenbouw.  Dit  betekent  dat  de  geschaalde  afstand  van 
de  C-zone  tweemaal  de  geschaalde  afstand  voor  huizenbouw  (SD  =  22,2)  bedraagt. 
De  afstand  van  de  C-zone  is  dus  afhankelijk  van  de  totale  ‘sympatische’  hoeveel¬ 
heid  explosief  volgens  de  volgende  relatie: 

D  =  44,4  (m) 

Hierbij  is  D  de  afstand  tot  de  donor  in  meter  en  Q  de  totale  hoeveelheid  opgeslagen 
explosief  in  kg  TNT. 

Deze  aaimame  wordt  later  via  de  output-file  op  waarheid  getoetst.  Hierbij  wordt 
geverifieerd  of  objecten  aan  de  rand  van  het  in  kaart  gebrachte  gebied  geen  signifi- 
cante  bijdrage  aan  het  individueel  risico  opleveren. 

Een  voorbeeld  van  een  input-file  voor  ‘RISKANAL’  is  gegeven  in  bijlage  B. 

Na  de  invoer  van  alle  donors  en  acceptors  geeft  het  programma  nog  een  aantal 
opties.  De  eerste  optie  die  volgt  op  de  invoer  van  de  objecten  is  de  mogelijkheid 
om  de  data  te  begrenzen  (‘rezoning’).  Als  een  data-file  een  complex  weergeeft  met 
een  groot  gebied  er  omheen  en  er  vindt  een  berekening  plaats  voor  dit  complex  met 
een  kleinere  hoeveelheid  opgeslagen  explosief,  dan  worden  de  omringende 
objecten  buiten  een  handmatig  te  specificeren  veiligheidsafstand  genegeerd. 
Vervolgens  bestaat  de  mogelijkheid  om  interactief  een  of  meerdere  data  van 
magazijnen  of  omringende  objecten  te  wijzigen.  Ook  kan  nog  een  nieuw  item 
worden  toegevoegd  of  een  bestaand  item  worden  verwijderd. 
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Indien  de  gebruiker  dit  wenst  bestaat  de  mogelijkheid  om  interactief  de  resultaten 
van  de  berekening  van  bet  interne  risico  te  wijzigen.  Dit  kan  noodzakelijk  zijn 
wanneer  een  specifiek  magazijn  of  een  specifieke  situatie  of  maatregel  niet  in  bet 
programma  is  verwerkt:  de  berekening  resulteert  bijvoorbeeld  in  sympatiscbe 
reacties  welke  in  realiteit  niet  ontstaan  of  vice  versa. 

Tenslotte  moeten  nog  enkele  vragen  beantwoordt  worden  om  de  output  naar  bet 
scberm  of  de  printer  te  controleren.  Er  kunnen  tevens  scbetsen  van  bet  complex, 
zowel  met  als  zonder  zijn  omgeving,  gemaakt  worden  op  bet  scberm  of  op  de 
printer. 

Voor  de  berekening  van  bet  interne  en  exteme  risico  moeten  twee  parameters  worden 
gedefinieerd.  De  eerste  bepaalt  welke  veiligbeidsafstand  volgens  AC/258  gebruikt 
moet  worden,  de  maximumafstand  of  korter.  In  een  aantal  gevallen  kan  een  kortere 
afstand  worden  gebruikt,  bijvoorbeeld  als  geen  primaire  of  erg  gevoelige  klasse  1.1 
explosieven  aanwezig  is  of  als  slecbts  kleine  kalibers  van  gevaren(sub)klasse  1.2 
aanwezig  zijn.  De  tweede  parameter  geeft  de  mogelijkbeid  om  de  berekening  uit  te 
voeren  met  of  zonder  de  gevolgen  van  fragmenten  en  brokstukken.  De  constructie  van 
een  magazijn  of  de  bescberming  van  een  omringend  object  kan  zodanig  zijn  dat  de 
effecten  van  fragmenten  verwaarloosd  kunnen  worden. 

Het  input-gedeelte  van  ‘RISKANAL’  wordt  afgesloten  met  enkele  berekeningen  en 
output-opties. 

•  De  vierde  coordinaat  van  elk  magazijn  wordt  berekend. 

•  De  begincondities  kunnen  worden  uitgeprint  op  scberm  of  printer  (bijlage  C 
geeft  een  voorbeeld). 

•  Scbetsen  van  de  beginsituatie  kunnen  worden  gemaakt  (opties:  alleen  de  maga- 
zijnen,  de  magazijnen  met  de  omgeving,  objecten  in  de  omgeving  weergegeven 
met  bet  symbool  X  of  met  een  pictogram,  dikke  of  dunne  lijnen,  A3  of  A4-for- 
maat)  (zie  bijlage  C). 

•  De  totale  boeveelbeid  betrokken  explosief,  onafbankelijk  van  de  gevaren(sub)- 
klasse,  wordt  berekend  met  bet  doel  objecten  die  buiten  de  veiligbeidscirkel, 
beborende  bij  deze  boeveelbeid  explosief  liggen  buiten  de  berekening  te  bou- 
den. 

•  Elk  magazijn  wordt  gekarakteriseerd  door  een  nummer,  afbankelijk  van  zijn 
mboud  (1  voor  gevaren(sub)klasse  1.1,2  voor  gevaren(sub)klasse  1.2,  3  voor 
gevaren(sub)klasse  1.3a,  4  voor  gevaren(sub)klasse  1.3b  en  5  voor  geva- 
ren(sub)klasse  1.4). 

‘RISKANAL’  is  nu  klaar  om  te  gaan  rekenen. 


2.2  Berekening  ‘intern  risico’ 

Een  complex  wordt  ‘intern  veilig’  genoemd  als  er  geen  sympatiscbe  reacties  optre- 
den  bij  een  explosie  van  een  willekeurig  magazijn.  Dit  deel  van  bet  programma 
ricbt  zicb  dus  op  de  berekening  van  mogelijke  sympatiscbe  reacties.  Explosieven- 
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opslagmagazijnen  dienen  op  een  zodanige  afstand  van  elkaar  en  van  andere  ge- 
bouwen  geplaatst  te  worden  dat  een  minimaal  haalbaar  risico  voor  mens,  dier  en 
bezit  resulteert.  In  de  ACy258-aanbevelingen  wordt  de  minimaal  toegestane  afstand 
tussen  een  donor  met  een  gegeven  hoeveelheid  explosieven  en  een  blootgesteld 
object  (acceptor)  gedefinieerd  als  de  ‘quantity  distance’  (QD).  De  ‘quantity  distan¬ 
ces’  zijn  gebaseerd  op  een  groot  aantal  tests  en  een  grondige  analyse  van  alle 
beschikbare  data  op  het  gebied  van  accidentele  explosies  in  verschillende  landen. 
Afhankelijk  van  het  type  blootgesteld  object  bestaan  verschillende  ‘quantity  dis¬ 
tances’.  Indien  het  blootgesteld  object  een  naastgelegen  magazijn  is,  wordt  de 
betreffende  QD  ook  wel  ‘Inter  Magazine  Distance’  genoemd. 

De  waarde  van  de  ‘quantity  distance’  is  echter  een  bron  van  onzekerheid  gezien  de 
grote  varieteit  in  explosies.  In  de  ACy258-richtlijn,  Part  I,  Annex  A,  Section  II 
worden  de  ‘quantity  distances’  gegeven  voor  de  diverse  situaties.  In  sommige 
gevallen  wordt  de  ‘quantity  distance’  gegeven  als  een  vaste  afstand,  in  andere 
gevallen  is  deze  afhankelijk  van  de  maximale  hoeveelheid  betrokken  explosief  in 
het  donormagazijn.  Aan  de  hand  van  deze  gegevens  uit  de  AC/258  berekent 
‘RISKANAL’  de  ‘quantity  distances’  (QD)  en  de  (kortste)  werkelijke  afstanden 
(DIST)  voor  elke  mogelijke  donor-acceptor  combinatie.  Voor  de  QD  van  igloos 
wordt  tevens  de  orientatie  van  de  magazijnen  in  rekening  gebracht.  Met  andere 
woorden:  ‘RISKANAL’  onderzoekt  of  andere  magazijnen  worden  geinitieerd 
(acceptormagazijnen)  door  66n  voor  66n  elk  magazijn  te  laten  detoneren  (donorma- 
gazijnen).  Alle  mogelijke  explosielocaties  met  zijn  mogelijke  sympathisch  reage- 
rende  magazijnen  worden  opgeslagen. 

Als  een  sympatische  reactie  optreedt,  kan  dit  acceptormagazijn  weer  sympatische 
reacties  veroorzaken  in  andere  magazijnen.  Er  worden  verschillende  rekenregels 
gebruikt  voor  de  verschillende  situaties,  hoofdzakelijk  in  overeenstemming  met 
AC/258.  In  tabel  2,  3  en  4  wordt  een  overzicht  gegeven  van  de  gebruikte  rekenre¬ 
gels  (hierin  is  QD  de  ‘quantity  distance’  en  CD  de  ‘conversion  distance’).  In  de 
AC/258  worden  verschillende  niveaus  aangegeven  voor  de  QD,  namelijk  Di  tot  en 
met  Di2.  De  ‘conversion  distance’  is  gedefinieerd  als  een  niveau  lager  dan  het  QD- 
niveau.  Als  bijvoorbeeld  voor  1.1  explosief  de  QD  gelijk  is  aan  D9  =  4.8Q^^3,  dan 
is  de  CD  gedefinieerd  als  Dg  =  3.6Q*^3.  De  term  ‘conversion  distance’  is  geintro- 
duceerd  om  te  kunnen  bepalen  of  explosieven  met  gevaren(sub)klasse  1.2  of  1.3  als 
gevaren(sub)klasse  1.1  behandeld  moeten  worden,  als  zij  te  dicht  bij  een  1.1 -bron 
zijn  opgeslagen  (volgens  de  regels  in  tabel  2,  3  en  4). 

Het  programma  telt  de  totale  hoeveelheid  explosieven  van  gevaren(sub)klasse  1.1 
voor  elke  mogelijke  donor  en  zijn  sympathisch  reagerende  acceptors  bij  elkaar  op. 
Het  centrum  van  de  explosie  (initiatiepunt)  van  de  1.1-detonatie  is  het  donormaga¬ 
zijn  zelf  als  dit  explosieven  van  een  gevaren(sub)klasse  1.1  bevat.  Is  dit  niet  het 
geval,  dan  bestaat  de  mogelijkheid  dat  een  of  meerdere  acceptormagazijnen  die  wel 
1.1 -explosieven  bevatten  sympathisch  zullen  reageren.  In  dat  geval  wordt  het 
acceptormagazijn  met  1.1 -explosieven  dat  het  dichtste  bij  de  donor  ligt  als  het 
centrum  van  de  explosie  gekozen. 
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Indien  alle  potentiele  explosielocaties  (donors)  op  deze  manier  zijn  berekend, 
wordt  een  label  met  de  verkregen  resultaten  op  bet  scherm  weergegeven.  Deze 
label  wordt  tevens  opgeslagen  in  een  data-file. 

Op  dit  punt  wordt  de  mogelijkheid  gecreeerd  om  handmatig  correcties  aan  te 
brengen  voor  situaties  welke  niet  kunnen  worden  beschreven  met  de  QD-theorie. 
Dit  kan  interactief  gedaan  worden,  waama  een  nieuwe  label  met  de  gewijzigde 
resultaten  wordt  gegeven. 

De  resultaten  van  de  berekening  van  bet  ‘interne  risico’  kunnen  als  volgt  worden 
samengevat. 

•  Voor  elke  potentiele  explosiebron  (donormagazijn)  worden  alle  sympatbiscb 
reagerende  magazijnen  (acceptors)  weergegeven. 

•  Voor  elk  donormagazijn  wordt  de  totale  boeveelbeid  betrokken  explosieven  van 
gevaren(sub)klasse  1.1  gegeven. 

•  Het  centrum  van  de  1 . 1  -explosie  wordt  gegeven. 

•  De  betrokken  boeveelbeid  explosieven  van  een  andere  gevaren(sub)klasse  en 
niet  reagerend  als  1.1  wordt  per  magazijn  gegeven. 


Tabel  2 :  Rekenregels  voor  een  1 . 1  -donor. 

Donor  D  gevaren- 

Acceptor  A  gevaren- 

Rekenregels 

(sub)klasse 

(sub)klasse 

1.1 

1.1 

A  binnen  1 .1  QD  van  D  geeft  een  sympatische  reactie  van  A 

1.1 

1.2/1 .3 

A  binnen  1 .1  QD  van  D  geeft  een  sympatische  reactie  van  A 

A  tevens  binnen  1 .1  CD  van  D  betekent  dat  de  inhoud  van  A  als 

1 .1-klasse  wordt  beschouwd. 

Tabel  3:  Rekenregels  voor  een  1.2-donor. 

Donor  D  gevaren- 

Acceptor  A  gevaren- 

Rekenregels 

(sub)klasse 

(sub)klasse 

1.2 

1.1 

A  binnen  1 .2  QD  van  D  geeft  een  sympatische  reactie  van  A 

A  tevens  binnen  1.1  CD  van  D  betekent  dat  de  inhoud  van  D  als 

1 .1-klasse  wordt  beschouwd 

1.2,  als  1.1 

1.1 

A  binnen  het  maximum  van  1 .1  QD  en  1 .2  QD  van  D  geeft  een 
sympatische  reactie  van  A 

1.2,  als  1.1 

1.2/1 .3 

A  binnen  het  maximum  van  1.1  en  1 .2  QD  van  D  geeft  een 
sympatische  reactie  van  A 

A  tevens  binnen  1 .1  CD  van  D  betekent  dat  de  inhoud  van  A  als 

1 .1-klasse  wordt  beschouwd. 

1.2 

1.2/1 .3 

A  binnen  1 .2  QD  van  D  geeft  een  sympatische  reactie  van  A 
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Tabel  4:  Rekenregels  voor  een  1.3 -donor. 


Donor  D  gevaren- 

Acceptor  A  gevaren- 

Rekenregels 

(sub)klasse 

(sub)klasse 

1.3 

1.1 

A  binnen  1.3  QD  van  D  geeft  een  sympatische  reactie  van  A 

A  tevens  binnen  1.1  CD  van  D  betekent  dat  50%  van  de  inhoud  van  D 

als  1.1-klasse  wordt  beschouwd  (gedeelte  verbrand  voordat  reactie 
plaatsvindt) 

1.3,  als  1.1 

1.1 

A  binnen  het  maximum  van  1.1  QD  en  1 .3  QD  van  D  geeft  een 
sympatische  reactie  van  A 

1.3,  als  1.1 

1.2/1 .3 

A  binnen  het  maximum  van  1.1  QD  en  1 .3  QD  van  D  geeft  een 
sympatische  reactie  van  A 

A  tevens  binnen  1.1  CD  van  D  betekent  dat  de  inhoud  van  A  als 

1 .1-klasse  wordt  beschouwd 

1.3 

1.2/1 .3 

A  binnen  1.3  QD  van  D  geeft  een  sympatische  reactie  van  A 

2.3  Berekening  ‘extern  risico’ 

Indien  alleen  de  berekening  van  het  ‘interne  risico’  gewenst  is,  kan  het  programma 
hier  worden  afgesloten.  Het  volgende  onderdeel  van  het  programma  berekent  het 
exteme  risico.  Hierbij  worden  de  fysische  explosie-effecten  gekwantificeerd,  de 
effecten  op  de  omgeving  berekend  en  de  gevolgen  voor  personen  in  de  omgeving 
(kans  op  letaliteit)  berekend. 

In  sommige  gevallen  reageren  groepen  van  magazijnen  sympathisch  en  op  een 
vergelijkbare  manier  en  geeft  elk  betrokken  magazijn  (donor)  dezelfde  exteme 
effecten.  Om  rekentijd  te  reduceren  kan  de  berekening  van  het  effect  in  dat  geval 
worden  weggelaten  door  deze  gelijk  te  stellen  aan  eerder  berekende  effecten. 

Als  onderdeel  van  de  berekening  zal  de  inhoud  van  elk  sympathisch  reagerend 
magazijn  worden  gedefinieerd.  De  tussenresultaten  van  de  berekeningen  kunnen  op 
het  scherm  of  op  een  file  worden  weergegeven. 

Om  de  fysische  effecten  te  kwantificeren  wordt  een  ‘loop’  uitgevoerd  voor  alle 
omiingende  objecten.  De  afstand  tussen  elk  object  en  het  reagerende  magazijn 
wordt  bepaald.  Voor  de  berekening  worden  wegen  opgedeeld  in  segmenten  van 
25  m,  waama  elk  segment  wordt  behandeld  als  een  apart  object,  gepositioneerd  op 
de  middencoordinaten  van  het  segment. 

Als  het  object  een  POL-faciliteit  betreft,  gelegen  birmen  25  of  60  m  van  het  maga¬ 
zijn,  dan  zal  een  bericht  hiervan  op  het  scherm  en  op  de  file  verschijnen.  Voor 
1.1-explosieven  dient  een  minimale  afstand  van  25  m  gehanteerd  te  worden.  Af- 
hankelijk  van  het  feit  of  de  drakontlasting  van  een  potentiele  explosielocatie  ge- 
richt  is  naar  de  POL-faciliteit  is  de  minimale  afstand  tot  deze  explosielocatie  voor 
1.3-explosieven  gelijk  aan  25  m  of  60  m  [2]. 

‘RISKANAL’  zal  vervolgens  de  explosie-effecten,  de  gevolgen  van  deze  effecten 
voor  de  omgeving  en  de  gevolgen  voor  personen  berekenen.  De  berekeningen  die 
hiema  worden  beschreven  worden  voor  elk  donormagazijn  gemaakt. 
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De  volgende  explosie-effecten  worden  in  ‘RISKANAL’  berekend: 

•  uitworp  van  fragmenten  en  brokstukken: 

-  aantal  fragmenten  per 

-  massaverdeling  en  eindsnelheid; 

•  ‘blast’  (druk  en  impuls): 

-  piekoverdruk  invallend; 

-  piekoverdruk  gereflecteerd; 

-  impuls  (invallend  en  gereflecteerd); 

-  positieve  faseduur, 

-  aankomsttijd  van  de  schokgolf; 

-  schokgolfsnelheid; 

•  hittestraling: 

-  intensiteit; 

-  duur  (blootstellingstijd). 

Voor  de  fragment-  en  brokstukdichtheid  worden  door  bet  TNO-PML  afgeleide 
formules  gebruikt,  gebaseerd  op  experimentele  data  uit  de  literatuur  [4,  5,  6].  Deze 
formules  werden  door  TNO  gevalideerd  door  berekende  resultaten  te  vergelijken 
met  meer  recente  data  [8, 9, 10],  waarbij  geen  significante  afwijkingen  werden 
geconstateerd.  De  fragment-  en  brokstukdichtheid  is  afhankelijk  van  het  type 
constructie  van  het  magazijn  en  van  de  hoeveelheid  opgeslagen  explosieven.  De 
overige  explosie-effecten  worden  berekend  met  formules  uit  de  ACy258-aanbeve- 
lingen.  Een  overzicht  van  alle  formules  is  in  bijlage  D  bijgevoegd. 

Na  de  berekening  van  de  fragmentdichtheid  worden  de  gevolgen  van  deze  frag¬ 
menten  en  brokstukken  bepaald.  Hiervoor  wordt  de  fragmentdichtheid  opgedeeld 
in  drie  gewichtsklassen: 

•  klasse  mi  zijn  alle  brokstukken  kleiner  dan  0,1  kg; 

•  klasse  m2  zijn  alle  brokstukken  tussen  0,1  kg  en  4,5  kg; 

•  klasse  m3  zijn  alle  brokstukken  groter  dan  4,5  kg. 

Voor  explosies  van  gevaren(sub)klasse  1.2  stoffen  wordt  aangenomen  dat  alle 
fragmenten  kleiner  zijn  deui  0,1  kg.  Als  functie  van  de  geschaalde  afstand  is  een 
verdeling  gemaakt  zoals  is  weergegeven  in  tabel  5. 

De  maximale  beginsnelheid  van  een  klein  fragment  (mi)  wordt  gesteld  op 
1500  m/s  (Gumey).  Voor  de  gemiddelde  beginsnelheid  van  een  klein  fragment  in 
de  richting  loodrecht  op  de  deur  van  een  igloo-type  magazijn  wordt  1200  m/s 
genomen.  Voor  een  niet  igloo-type  magazijn  wordt  de  gemiddelde  beginsnelheid 
aan  de  deurzijde  gesteld  op  1000  m/s.  Aan  de  zijkant  en  achterzijde  van  het  maga¬ 
zijn  worden  deze  snelheden  gesteld  op  respectievelijk  500  m/s  en  300  m/s.  Voor 
fragmenten  en  brokstukken  uit  klasse  m2  wordt  een  beginsnelheid  van  300  m/s 
gehanteerd  aan  de  ingang  van  een  igloo-type  magazijn.  Voor  de  overige  situaties 
wordt  een  beginsnelheid  van  200  m/s  aangenomen.  Voor  de  zware  fragmenten 
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(m3)  tenslotte,  wordt  een  beginsnelheid  aangenomen  van  200  m/s  aan  de  ingang 
van  een  igloo-type  en  een  snelheid  van  100  m/s  in  de  overige  situaties. 


Tabel  5:  Massaverdeling  van  fragmenten  en  brokstukken  alsfunctie  van  de  geschaal- 

de  afstand. 


Geschaalde  afstand 

(m/kg^^3) 

mi 

Massaverdeling  (%) 

1712 

m3 

SD<1 

75 

20 

5 

1  <SD<2 

70 

26 

4 

2<SD<3 

65 

31 

4 

3<SD<4 

60 

36 

4 

4<SD<5 

55 

41 

4 

5<SD<6 

50 

47 

3 

6  <  SD  <  7 

45 

52 

3 

7<SD<8 

40 

57 

3 

8<SD 

35 

62 

3 

De  snelheid  op  een  bepaalde  afstand  x  van  de  donor  (objectafstand)  wordt  bere- 
kend  met  de  formule: 

Vx  =  Vo 

Met: 

Vx  =  snelheid  op  afstand  x  (m/s) 

Vq  =  beginsnelheid  (m/s) 

L  =  een  parameter  gedefinieerd  door  L  =  2{k^  •  m)l^3  /  .p)  ^let: 

k  =  vormfactor  (kg/m^) 

m  =  firagmentmassa  (kg) 

Cd  =  dragcoefficient 
p  =  dichtheid  van  het  medium  (kg/m^) 

Als  gelijkvormige  fragmenten  met  een  vormfactor  k  =  2,61-103  kg/m^ 

(2,61  g/cm3),  een  dragcoefficient  Cd  =  1,28  en  de  dichtheid  van  lucht  worden 
aangenomen,  levert  dit: 

L  =  247  (m) 

voor  gietstalen  projectielen  en  ffagmentatiebommen  in  lucht  bij  standaardcondities 
en  verwaarlozing  van  de  zwaartekracht  (AC/258,  part  II). 

Voor  penetratie  in  beton  wordt  vervolgens  de  maximale  penetratiediepte  berekend 
met  de  volgende  formule: 

X  =  1, 0764-106  .  ^0.4  .  y^,8  +  0,04752  ■  rrfi-33  (m) 


Als  een  gemiddelde  en  constante  muurdikte  d  van  0,1  m  wordt  aangenomen  kan  de 
restsnelheid  na  penetratie  berekend  worden  met: 
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Vr=V ^^1-075^  (mis) 

Voor  1,5  mm  staal  (auto’s)  en  voor  de  menselijke  huid  wordt  de  snelheid  berekend 
waarop  fragmenten  van  respectievelijk  0,1, 1,0  en  4,5  kg  een  kans  hebben  van  50% 
dat  penetratie  plaatsvindt  (V50  ).  Dit  is  de  zogenaamde  ‘ballistic  limit  velocity’.  De 
formules  die  hiervoor  gebruikt  worden  zijn  gegeven  in  bijlage  E.  Met  de  snelheden 
Vx  en  V50  kan  de  restsnelheid  na  perforatie  van  de  beschermende  laag  (huid,  staal), 
Vr  berekend  worden: 

Vr=l(Vl-vio) 

Deze  restsnelheid  is  dan  bepalend  voor  de  berekening  van  de  kans  op  letaliteit. 

Als  alle  fysische  parameters  zijn  bepaald  worden  de  effecten  op  mensen  voorspeld 
door  middel  van  probit-functies.  Deze  probit-functies  zijn  door  TNO-PML  afgeleid 
van  experimentele  data  uit  de  literatuur  aangevuld  met  data  verkregen  uit  eigen 
experimenten.  De  probit-functies  zijn  beschreven  in  het  ‘Groene  Boek’  [3]  waarin 
kwetsbaarheidsmodellen  voor  de  effecten  op  mensen  zijn  gegeven.  De  algemene 
vorm  van  deze  probit-functies  is: 

Pr  =  a  +  b  ln  S 

Hierbij  zijn  a  en  b  constanten  en  is  S  een  variabele.  Afhankelijk  van  de  waarde  van 
S,  geeft  de  grootte  van  Pr  de  kans  op  een  bepaalde  gebeurtenis  (voor  dit  geval 
hoofdzakelijk  ae  kans  op  letaliteit).  De  kans  (in  %)  behorende  bij  de  Pr-waarden  is 
weergegeven  in  tabel  6. 

Tabel  6:  Relatie  tussen  kansen  en  probits. 


% 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

0,00 

2,67 

2.95 

3,12 

3,25 

3,36 

3,45 

3,52 

3,59 

3.66 

10 

3,72 

3.77 

3,82 

3,897 

3,92 

3,96 

4.01 

4,05 

4,08 

4,12 

20 

4,16 

4,19 

4,23 

4,26 

4,29 

4,33 

4,36 

4,39 

4,42 

4,45 

30 

4.48 

4,50 

4,53 

4,56 

4,59 

4,61 

4,64 

4,67 

4,69 

4.72 

40 

4,75 

4,77 

4,80 

4,82 

4,85 

4,87 

4,90 

4,92 

4,95 

4,97 

50 

5.00 

5,03 

5,05 

5,08 

5.10 

5,13 

5,15 

5,18 

5.20 

5,23 

60 

5,25 

5,28 

5,31 

5.33 

5,36 

5,39 

5,41 

5.44 

5.47 

5,50 

70 

5,52 

5.55 

5,58 

5,61 

5,64 

5,67 

5.71 

5,74 

5,77 

5,81 

80 

5,84 

5,88 

5,92 

5,95 

5,99 

6,04 

6,08 

6,13 

6,18 

6,23 

90 

6,28 

6,34 

6,41 

6,48 

6.55 

6,64 

6,75 

6,88 

7,05 

7,33 

Als  Pr  kleiner  dan  0  is  wordt  aangenomen  dat  de  kans  gelijk  is  aan  nul  en  als  Pr 
groter  is  dan  7,33  wordt  de  kans  op  100%  gesteld.  Voor  een  Pr-waarde  van  4,88 
bijvoorbeeld,  wordt  de  bijbehorende  kans  van  46%  genomen.  Txt  worden  de  vol- 
gende  effecten  op  mensen  berekend: 

•  de  effecten  als  gevolg  van  fragmenten  en  brokstukken: 

-  fragmenten  <0,1  kg:  criterium  huidperforatie; 

-  fragmenten  tussen  0,1  kg  en  4,5  kg:  criterium  kinetische  energie; 
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-  fragmenten  >  4,5  kg:  criterium  schedelbasisfractuur; 

•  de  effecten  als  gevolg  van  ‘blast’  (druk  en  impuls): 

-  trommelvliesbreuk  (=  niet  letaal); 

-  letale  longschade; 

-  letaal  letsel  als  gevolg  van  schade  aan  of  instorting  van  gebouwen; 

-  letaal  letsel  als  gevolg  van  ruitbreuk; 

-  letaal  letsel  door  meesleureffect  als  gevolg  van  ‘blast’; 

•  de  effecten  als  gevolg  van  hittestraling: 

-  2^  graads  verbranding; 

-  3®  graads  verbranding. 

Hierbij  worden  de  volgende  aannames  gedaan. 

•  Voor  de  effecten  van  fragmenten  wordt  voor  elke  massaklasse  berekend  wat  de 
kans  is  op  letaliteit  als  een  mens  door  een  dergelijk  fragment  wordt  geraakt.  Het 
totale  oppervlak  van  een  mens  wat  bij  inslag  tot  letaliteit  zou  leiden  wordt  ge- 
schat  op  0,58  m^  voor  rechtopstaande  personen.  Dit  oppervlak  wordt  gebruikt 
om  de  kans  op  raken  te  berekenen. 

•  Voor  de  effecten  van  ruitbreuk  wordt  de  kans  op  letaliteit  opgebouwd  uit  de 
kans  op  ruitbreuk,  de  kans  dat  ruitbreuk  leidt  tot  letaliteit  bij  personen  binnen 
1,75  m  achter  de  ruit  en  de  kans  dat  een  persoon  zich  in  de  ruimte  binnen 

1,75  m  achter  de  ruit  bevindt.  Voor  de  kans  op  ruitbreuk  worden  twee  verschil- 
lende  probit-functies  gebruikt,  namelijk  voor  oude  huizen  (gebouwd  voor  1975) 
en  nieuwe  huizen  (gebouwd  na  1975).  Voor  de  kans  dat  een  persoon  zich  bin¬ 
nen  1,75  m  achter  een  ruit  bevindt  wordt  5%  gehanteerd.  Deze  waarde  is  relatief 
laag  voor  kantoren  tijdens  kantooruren,  maar  relatief  hoog  voor  huizen. 

Voor  een  overzicht  van  alle  probit-functies  die  gebruikt  worden  in  ‘RISKANAL’ 
wordt  verwezen  naar  het  ‘Groene  Boek’  [3]  of  van  Amelsfort  [1]. 

Aangezien  de  formules  die  door  ‘RISKANAL’  gebruikt  worden  slechts  geldig  zijn 
voor  relatief  grote  hoeveelheden  explosieven  levert  dit  conservatieve  waarden  op 
voor  de  kans  op  letaliteit  indien  slechts  kleine  hoeveelheden  explosieven  zijn 
opgeslagen.  Nadat  alle  effecten  zijn  berekend,  wordt  daarom  door  het  programma 
een  correctie  gemaakt  voor  de  opslag  van  kleine  hoeveelheden  explosieven  (kleiner 
dan  500  kg).  In  dat  geval  worden  de  aangepaste  veiligheidsafstanden  gehanteerd 
als  gegeven  in  [14]  en  in  tabel  7. 

Als  een  omringend  object  geplaatst  is  tussen  de  afstanden  voor  ‘geen  gevaar’  en 
‘lichte  bouwkundige  schade’,  wordt  de  kans  op  letaliteit  resulterend  uit  de  bereke- 
ningen  gedeeld  door  tien.  Voor  objecten  buiten  de  afstand  voor  ‘geen  gevaar’ 
wordt  de  kans  op  letaliteit  op  nul  gesteld. 
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Tabel  7:  Aangepaste  veiligheidsafstanden  voor  opslag  van  kleine  hoeveelheden  explosieven. 


Hoeveelheid 

(kg) 

Inhabited  building  distances  (m) 
gevaren(sub)klassen  1.1  en  1.2 
geen  gevaar  lichte  bouwkundige  schade 

gevaren(sub)klasse  1 .3 

<50 

75 

45 

15 

51  - 100 

105 

65 

15 

101 -150 

125 

75 

25 

151 -200 

135 

75 

25 

201  -  250 

145 

85 

25 

251  -  300 

155 

85 

30 

301  -  350 

165 

95 

30 

351  -  400 

175 

95 

30 

401  -  450 

175 

105 

35 

451  -  500 

185 

105 

35 

Nu  zijn  alle  kansen  op  letaliteit  voor  alle  objecten  en  alle  bronnen  voor  een  bepaal- 
de  donor  berekend.  Er  kan  een  statistische  optelsom  van  alle  kansen  voor  alle 
objecten  worden  gemaakt.  Hierbij  worden  de  volgende  aannames  gedaan. 

•  Het  aantal  gewonden  als  gevolg  van  instorting  van  gebouwen  wordt  gesteld  op 
twee  keer  het  aantal  letaliteiten. 

•  Het  aantal  gewonden  als  gevolg  van  fragmenten  wordt  gesteld  op  twee  keer  het 
aantal  letaliteiten. 

•  Het  aanteil  gewonden  als  gevolg  van  het  meesleureffect  als  gevolg  van  blast 
wordt  gesteld  op  twee  keer  het  aantal  letaliteiten. 

De  kans  op  letaliteit  wordt  voor  elk  object  en  voor  de  bron  op  het  scherm  weerge- 
geven  en  naar  een  file  geschreven.  Nu  kunnen  de  effecten  van  alle  bronnen  van 
deze  donor  voor  alle  objecten  berekend  worden.  Voor  wegen  wordt  een  gemiddel- 
de  waarde  voor  de  kans  op  letaliteit  genomen  voor  auto’s  en  fietsers/voetgangers, 
waarbij  een  gemiddelde  snelheid  van  60  km/uur  voor  auto’s  en  15  km/uur  voor 
fietsers/voetgangers  wordt  aangenomen. 

De  kans  op  letaliteit  voor  een  bepaald  object  is  de  kans  op  letaliteit  maal  het  aantal 
units  maal  het  aantal  personen  per  unit  maal  de  gemiddelde  bezetting  van  het 
object.  Voor  wegen  wordt  per  segment  van  25  m  de  gemiddelde  bezetting  bepaald, 
welke  afhankelijk  is  van  de  verkeersintensiteit  en  de  snelheid.  Voor  elk  segment 
wordt  de  kans  op  letaliteit  vermenigvuldigd  met  de  verkeersbezetting,  waama  deze 
getallen  voor  alle  segmenten  bij  elkaar  worden  opgeteld.  Bijlage  F  geeft  een  voor- 
beeld  van  een  standaard  gedetailleerde  output  zoals  deze  door  het  programma 
wordt  gegeven. 


2.4  Output 

In  de  output  file  worden  de  effecten  van  alle  mogelijke  donors  van  het  complex 
weergegeven  en  volgt  een  overzicht  van  de  risico’s  voor  het  complex.  Voor  elke 
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donor  wordt  weergegeven  welk  type  donor  het  betreft,  wat  de  massa  en  geva- 
ren(sub)klasse  van  de  inhoud  is  en  voor  elk  object  de  effecten  als  gevolg  van  een 
initiatie  in  deze  donor.  Per  object  worden  de  volgende  effecten  aangegeven: 

•  piekoverdruk; 

•  impuls; 

•  fragmentdichtheid; 

•  positieve  fase-duur; 

•  stralingsintensiteit; 

•  blootstellingstijd; 

•  kans  op  trommelvliesbreuk; 

•  kans  op  longschade; 

•  kans  op  letaliteit  door  fragmenten; 

•  kans  op  letaliteit  door  fragmenten  voor  fietsers/voetgangers; 

•  kans  op  letaliteit  door  instorting  van  huizen; 

•  kans  op  letaliteit  door  ruitbreuk; 

•  kans  op  letaliteit  door  meesleureffect  als  gevolg  van  blast; 

•  kans  op  letaliteit  door  hittestraling  (3®  graads  verbranding); 

•  kans  op  2®  graads  verbranding. 

Per  donor  wordt  de  output  afgesloten  met  een  tabel  met  daarin  per  object  de  kans 
op  letaliteit  en  de  kans  op  verwonding.  Voor  wegen  worden  de  kansen  voor  auto- 
mobilisten  en  voor  fietsers/voetgangers  apart  aangegeven. 

Vervolgens  worden  voor  elk  object  de  kansen  op  letaliteit  voor  elke  donor  bij 
elkaar  opgeteld  en  gedeeld  door  het  aantal  donors.  Dit  geeft  de  gemiddelde  kans  op 
letaliteit  per  object  voor  dit  complex  en  wordt  in  een  standaardtabel  in  de  output- 
file  weergegeven.  Door  deze  gemiddelde  kans  op  letaliteit  te  vermenigvuldigen 
met  de  kems  op  een  explosie  wordt  het  gemiddelde  individueel  risico  (IR)  per 
object  gevonden.  Deze  is  echter  niet  in  de  tabel  opgenomen.  In  dezelfde  tabel 
wordt  wel  aangegeven  welke  bron  de  grootste  kans  op  letaliteit  geeft  en  hoe  groot 
die  kans  voor  deze  bron  is  (bijlage  G). 

Uit  de  gemiddelde  kans  op  letaliteit  per  object  en  het  aantal  personen  verbonden 
aan  dit  object  kan  een  gemiddeld  aantal  letaliteiten  per  object  worden  berekend. 
Deze  gemiddelde  letaliteiten  per  object  bij  elkaar  opgeteld  geeft  het  gemiddeld 
aantal  te  verwachten  letaliteiten  voor  het  complex.  Dit  getal  wordt  in  de  output-file 
weergegeven.  Bovendien  wordt  het  gemiddelde  aantal  te  verwachten  zware  ver- 
wondingen  weergegeven. 

Tenslotte  wordt  het  groepsrisico  bepaald  door  het  gemiddelde  aantal  te  verwachten 
letaliteiten  voor  het  complex  te  vermenigvuldigen  met  de  kans  op  een  explosie.  De 
kans  op  een  explosie  is  gelijk  aan  de  kans  op  een  initiatie  in  een  bepaald  magazijn 
(=  10'^)  maal  het  aantal  magazijnen.  Ook  het  groepsrisico  wordt  in  de  output-file 
weergegeven. 
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Het  prograrmna  ‘RISKANAL’  geeft  nu  nog  de  mogelijkheid  om  een  plaatje  te 
genereren  op  het  scherm  of  op  een  plotter  van  alle  objecten  waarbij  letaliteit  of 
verwonding  te  verwachten  valt. 

Het  programma  wordt  vervolgens  afgesloten  met  een  aantal  vragen. 

•  Wilt  u  een  nieuwe  data-file  maken? 

•  Wilt  u  een  nieuwe  ‘run’  uitvoeren? 

•  Indien  ja,  wilt  u  de  input-file  opnieuw  inlezen? 

Indien  op  al  deze  vragen  negatief  wordt  geantwoord  kan  het  programma  verlaten 
worden. 
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3  Conclusies 


In  dit  rapport  is  een  gedetailleerd  overzicht  gegeven  van  het  interactieve  program- 
ma  RISKANAL.  Dit  programma  wordt  gebruikt  om  het  interne  en  exteme  risico 
voor  een  munitiemagazijnencomplex  te  berekenen.  Het  programma  is  gebaseerd  op 
een  combinatie  van  de  NATO  AC/258-aanbevelingen  en  de  meest  recente  letsel- 
criteria  zoals  beschreven  in  het  Nederlandse  ‘Groene  Boek’.  De  ACy258- 
aanbevelingen  zijn  gebruikt  voor  de  berekening  van  het  intern  risico  van  het  com¬ 
plex  en  voor  de  berekening  van  de  explosie-effecten.  Voor  de  berekening  van  de 
fragment-  en  brokstukdichtheid  wordt  gebruikgemaakt  van  door  het  TNO  Prins 
Maurits  Laboratorium  afgeleide  formules,  gebaseerd  op  de  ESKIMO  en  Yuma- 
tests.  De  formules  uit  het  ‘Groene  Boek’  worden  gebruikt  voor  de  berekening  van 
het  extern  risico.  Hoewel  de  bruikbaarheid  van  enkele  criteria  gelimiteerd  is,  lijken 
de  eindresultaten  van  de  berekeningen  erg  realistisch. 

Aangezien  de  kennis  op  het  gebied  van  risicoanalyse  aan  ontwikkeling  onderhevig 
is,  is  het  programma  zo  opgezet  dat  het  in  de  toekomst  mogelijk  is  om  gebruikte 
relaties  eenvoudig  te  vervangen  door  meer  recente  relaties. 

Door  het  interactieve  karakter  van  het  programma  kunnen  de  input-gegevens 
eenvoudig  worden  veranderd  en  kunnen  de  effecten  van  deze  veranderingen  snel 
zichtbaar  worden  gemaakt.  Hierdoor  kan  op  een  relatief  eenvoudige  en  snelle 
manier  een  ontwerp  voor  een  optimale  lay-out  worden  gemaakt. 

De  absolute  getallen  dienen  met  enige  voorzichtigheid  te  worden  geinterpreteerd. 
Doordat  wereldwijd  slechts  beperkte  experimentele  data  beschikbaar  zijn,  zijn  veel 
van  de  formules  gebruikt  voor  risicoanalyse  slechts  beperkt  gevalideerd.  De  on- 
nauwkeurigheid  van  de  resultaten  kan  door  de  gebruikte  formules  relatief  groot 
zijn.  Om  dit  te  verbeteren  is  experimenteel  onderzoek  noodzakelijk. 

Echter,  de  verkregen  kwantitatieve  resultaten  zijn  een  goed  hulpmiddel  om  waar- 
den  voor  het  extern  risico  relatief  met  elkaar  te  vergelijken. 
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Bijlage  A  Effect  als  fiinctie  van  de  afstand  voor  verschil 
lende  hoeveelheden  explosieve  stoffen  (1.1) 


A.1 
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PPOb.  Df  injury  (>)  Ppob.  of  injury  (,) 


Figuur  A.2:  Effect  (kans  op  verwonding)  alsfunctie  van  de  afstand  voor  drie  hoeveelhe- 
den  explosieven  van  gevaren(sub)klasse  1.1,  inclusiefhet  effect  vanfrag- 
menten  en  brokstukken. 
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EffGCl  Of  1.1  eXpiDSiUEE. 


Effect  of  1.1  exploBides. 


Figuur  A.3:  Effect  (kans  op  letaliteit)  alsfunctie  van  de  afstand  voor  drie  hoeveelheden 
explosieven  van  gevaren(sub)klasse  1.1,  zonder  het  effect  vanfragmenten  en 
brokstukken. 
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Bijiage  B  Data-file  format  met  een  voorbeeld 

Card  1: 

IX 

number  of  magazines 

14 

number  of  surrounding  objects 

14 

number  of  object  coordinates 

14 

- 

IX 

average  number  of  people  in  houses 

F5.2 

- 

IX 

complex  name 

A50 

For  each  magazine: 

Card  2: 

IX 

magazine  code 

A2 

- 

2X 

magazine  name 

A6 

- 

2X 

explosive  contents  (kg) 

FIO.O 

- 

2X 

hazard  division 

A4 

Card  3: 

IX 

3  comer  coordinates 

6F12.5 

For  each  surrounding  object: 

Card  4: 

IX 

object  name 

A6 

- 

IX 

object  code 

A2 

- 

IX 

first  coordinate  number 

14 

last  coordinate  number 

14 

number  of  houses/road  intensity 

16 

-/cyclist  euid  pedestrian  intensity 

14 

- 

IX 

object  centre  or  first  road  coordinates  2F12.5 

IX 

6F12.5 


B.1 


Card  5; 


next  road  coordinates 
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Voorbeeld  van  een  input-file  voor  RISKANAL. 


5 

14 

20  3.25 

EXAMPLE 

13 

IGLOOl 

15000. 

1.1 

112.5 

125.0 

122.5 

125.0 

122.5 

150.0 

OB 

EARTH2 

50000. 

1.3a 

250.0 

100.0 

250.0 

90.0 

275.0 

90.0 

13 

IGL003 

15000. 

1.1 

250.0 

37.5 

250.0 

27.5 

275.0 

27.5 

OP 

OPEN4 

20000. 

1.2 

250.0 

- 

25.0 

250.0 

-35.0 

275.0 

-35.0 

13 

IGL005 

15000. 

1.1 

175.0 

-110.0 

175.0 

-120.0 

200.0 

-120.0 

HSl 

HS 

1 

1 

5 

0 

75.0 

600.0 

HS2 

HS 

2 

2 

8 

0 

150.0 

690.0 

HS3 

HS 

3 

3 

5 

0 

230.0 

600.0 

CM4 

HU 

4 

4 

90 

0 

550.0 

200.0 

CM5 

HU 

5 

5 

90 

0 

800.0 

700.0 

HS6 

HS 

6 

6 

5 

0 

420.0 

250.0 

PK7 

HU 

7 

7 

5 

0 

950.0 

140.0 

HS8 

HS 

8 

8 

5 

0 

420.0 

50.0 

HS9 

HS 

9 

9 

2 

0 

190.0 

-245.0 

HSIO 

HS 

10 

10 

1 

0 

10.0 

-245.0 

HSll 

HS 

11 

11 

1 

0 

-150.0 

50.0 

RDl 

RD 

12 

13 

1000 

50 

-450.0 

650.0 

1000.0 

650.0 

RD2 

RD 

14 

18 

100 

10 

-50.0 

650.0 

-50.0 

-200.0 

350.0  -200.0  450.0 

-100.0 

1000.0  -100.0 
RD3  RD  19  20 

395.0  700.0 


B.2 


10  100 


395.0 


-160.0 
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B.3 


Voorbeeld  van  hoe  het  programma  de  input-file  weergeeft. 

Berekeningen  met  maximale  ‘quantity  distances’. 
Berekeningen  met  gevaar  van  fragmenten. 


The  following  data  concern  the  storage  sites: 


IGLOOl 

Type  13 

High  explosive  contents 

15000.  kg 

Hazard  Division  1.1 

112.50 

125.00 

122.50  125.00  122.50 

150.00 

112.50 

150.00 

EARTH2 

Type 

OB 

High  explosive  contents 

50000.  kg 

Hazard  Division  1.3a 

250.00 

100.00 

250.00  90.00  275.00 

90.00 

275.00 

100.00 

IGL003 

Type  13 

High  explosive  contents 

15000.  kg 

Hazard  Division  1.1 

250.00 

37.50 

250.00  27.50  275.00 

27.50 

275.00 

37.50 

OPEN4 

Type  OP 

High  explosive  contents 

20000.kg 

Hazard  Division  1.2 

250.00 

-25.00 

250.00  -35.00  275.00 

-35.00 

275.00 

-25.00 

IGL005 

Type  13 

High  explosive  contents 

15000.  kg 

Hazard  Division  1.1 

175.00 

-110.00 

175.00  -120.00  200.00 

-120.00 

200.00 

-110.00 

The  following  data  concern  the  surroundings: 


Average  number  of  people  in  houses  (units)  is  3.25 


HSl 

Type  HS 

Number  of  houses 

5 

75.00 

600.00 

HS2 

Type  HS 

Number  of  houses 

8 

150.00 

690.00 

HS3 

Type  HS 

Number  of  houses 

5 

230.00 

600.00 

CM4 

Type  HU 

Number  of  houses 

90 

550.00 

200.00 

CM5 

Type  HU 

Number  of  houses 

90 

800.00 

700.00 

HS6 

Type  HS 

Number  of  houses 

5 

420.00 

250.00 

PK7 

Type  HU 

Number  of  houses 

5 

950.00 

140.00 

HS8 

Type  HS 

Number  of  houses 

5 

420.00 

50.00 

HS9 

Type  HS 

Number  of  houses 

2 

190.00 

-245.00 

HSIO 

Type  HS 

Number  of  houses 

1 

10.00 

-245.00 

HSll 

Type  HS 

Number  of  houses 

1 

-150.00 

50.00 
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B.4 


RDl  Type  RD 

-450.00 
1000.00 

RD2  Type  RD 

-50.00 
-50.00 
350.00 
450.00 
1000.00 

RD3  Type  RD 

395.00 
395.00 


24  Hours  intensity  1000 
650.00 
650.00 

24  Hours  intensity  100 

650.00 
-200.00 
-200.00 
-100.00 
-100.00 

24  Hours  intensity  10 

-160.00 
700.00 


Number  of  cyclists  50 


Number  of  cyclists  10 


Number  of  cyclists  100 
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Bijlage  C  Weergave  van  begincondities 


I 

EXflmPLE 

IGL001 


^EARTHS 
^IGL003 
OPEn4 

_ ,  IGL005 

SCALE  1  :  2700 

Figuur  C.l :  Voorbeeld  van  een  grafische  weergave  van  een  munitiemagazijnencomplex. 


SCALE  1  :  10500 


400  m 


Figuur  C.2:  Voorbeeld  van  de  grafische  weergave  van  het  munitiemagazijnencomplex  met 
zijn  omgeving. 
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D.1 


Bijiage  D  Formules  voor  het  kwantificeren  van  de  fysi- 
sche  effecten 


D.l  Fragmentdichtheid  (1.1  en  1.2-bron) 

De  formules  die  door  ‘RISKANAL’  worden  gebruikt  zijn  afgeleid  van  de  gefitte 
experimentele  data  verkregen  tijdens  de  ESKIMO  1  (M107-granaten,  10.900  kg, 
Weals,  1973),  ESKIMO  2  (M117-bommen,  10.900  kg.  Weals,  1974)  en  Yuma 
(M117-bommen,  2970  kg,  Feinstein,  1972)  testprogramma’s  [4,  5,  6]. 

De  berekende  resultaten  voor  de  fragmentdichtheid  op  een  bepaalde  afstand,  ver¬ 
kregen  uit  deze  data,  is  niet  in  tegenspraak  met  meer  recente  experimentele  studies 
zoals  de  ‘joint  Australian/UK  stack  fragmentation  trials:  phase  1’  in  1982,  phase  2 
in  1985  en  phase  3  in  1987  [8,  9, 10]. 

De  berekeningen  zijn  uitgevoerd  met  de  totale  hoeveelheid  1.1-explosieven.  De 
orientatie  van  de  igloos  ten  opzichte  van  het  object  wordt  meegenomen  in  de 
berekening. 

Gebruikte  symbolen: 

X  =  afstand  in  m; 

Q  =  hoeveelheid  explosieven  betrokken  bij  explosie  in  kg; 
q  =  fragmentdichtheid  in  l/m^. 

Extrapoleren  naar  waarden  van  q  boven  0,5  per  vierkante  meter  is  riskant.  Indien 
dit  zich  voordoet  wordt  een  kwadratische  toename  van  q  met  afnemende  afstand  x 
aangenomen. 

Fragmentdichtheid  voor  niet  igloo-type  magazijnen  met  minder  dan  3000  kg: 
log  q  =  9,7591  -  5,1203  log  (x/Q0.0104) 
als  q  >  0,5:  Z  =  (11, 484- IO^.QO, 541)0, 1953 

q  =  0,5  (Z/x)2 

Fragmentdichtheid  voor  niet  igloo-type  magazijnen  met  meer  dan  3000  kg: 
log  q  =  6,00252  -  5,08955  log  (x/QO’306) 
als  q  >  0,5:  Z  =  (2,012-106.Qi.557)0,l9648 

q  =  0,5  (Z/x)2 

Fragmentdichtheid  aan  de  zijkant  van  een  igloo-type  magazijn: 
log  q  =  5,64185  -  5,3782  log  (x/Q0.306) 
voor  een  1.2-donor  geldt:  q  =  0 
als  q  >  0,5:  Z  =  (8,7676-105.Ql.6457)0, 18594 

q  =  0,5  (Z/x)2 
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Fragmentdichtheid  aan  de  achterzijde  van  een  igloo-type  magazijn: 
log  q  =  4,10427  -  4,43798  log  (x/Q0.306) 
voor  een  1.2-donor  geldt:  q  =  0 
als  q  >  0,5:  Z  =  (2,5428- 104.Qi.358)0, 22533 

q  =  0,5  (Z/x)2 

Fragmentdichtheid  aan  de  voorkant  van  een  igloo-type  magazijn: 
log  q  =  6,00252  -  5,08955  log  (x/Q0.306) 
als  1.1 -donor  en  q  >  0,5: 

Z  =  (2,012-106.Q1.557)0, 19648 

q  =  0,5  (Z/x)2 

Indien  een  magazijn  is  voorzien  van  een  aarden  omwalling  geldt  in  de  richting 
waarin  de  wal  staat  een  reductiefactor  0,5:  ofwel  q’=  q/2. 

Voor  1.2-opslag  wordt  uitgegaan  van  de  volgende  formules: 

als  Q  <  100  kg:  log  q  =  2,004722  -  0,030369  x  -i-  0,00006  x^  -  4,41 1449- 10-8  x3 

als  100  <  Q  <  500  kg: 

log  q  =  2,159189  -  0,027044  x  +  0,00005  x2  -  3,41991-10-8  x3 
als  500  <  Q  <  2000  kg: 

log  q  =  2,64663  -  0,027371  x  -i-  0,000056  x2  -  4,05308 1  10-8  x3 
als  Q  >  2000  kg: 

log  q  =  3,066205  -  0,025275  x  -i-  0,000047  x2  -  3,23537410-8  x3 

Voor  de  zijkant  en  achterkant  van  een  igloo  wordt  aangenomen  dat  q  =  0. 

Indien  een  magazijn  is  voorzien  van  een  aarden  omwalling  geldt  in  de  richting 
waarin  de  wal  staat  een  reductiefactor  1/5:  ofwel  q’=  q/5. 


D.2  Overige  fysische  effecten 

Alle  hiema  genoemde  effecten  worden  beschreven  met  de  gefitte  formules  als 
gegeven  in  de  AC/258  ‘Manual  on  NATO  safety  principles  for  the  storage  of 
ammunition  and  explosives’,  1977,  Part  II,  724.c,  april  1987. 

Gebruikte  symbolen: 

Ps  =  invallende  druk  in  kPa 
Pr  =  gereflecteerde  druk  in  kPa 
Is  =  invallende  impuls  in  kPas 
If  =  gereflecteerde  impuls  in  kPas 
Ta  =  schokgolfaankomsttijd  in  s 


D.2 
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Vs  =  snelheid  van  de  schokgolf  in  m/s 
Tp  =  duur  van  de  positieve  fase  in  s 
Rc  =  straal  van  de  vuurbal  in  m 
Rg  =  straal  van  de  vuurbal  op  grondniveau  in  m 
W  =  hittestralingsintensiteit  in  kW/m^ 
t  =  blootstellingstijd  in  s 
Q  =  betrokken  hoeveelheid  explosieven  in  kg 
X  =  afstand  in  m 

sc  =  scaled  distance  in  m/kg^^^r  sc  =  x/Q^^3 

Druk  (1.1-bronnen) 

Voor  open  opslag  en  lichte  constructies  worden  relaties  volgens  de  ACy258- 
richtlijnen  gebruikt.  Voor  igloos  en  zware  constructie  donormagazijnen  wordt 
gebruikgemaakt  van  NATO  US(ST)IWP/103-95:  ‘Airblast  from  earth-covered 
magazines:  incorporation  of  additional  data  and  attempts  to  reduce  the  scatter’. 

Open  opslag  en  lichte  constructies: 
sc  <  0,0674:  Ps  =  10^ 

so  40:  Ps  =  0 

0,0674  <  sc  <  40: 

U=  -0,214362789151  +  1,35034249993  log(sc) 

Y  =  +  2,78076916577  -  1,6958988741  U  -  0,154159376846 

+  0,514060730593  +  0,0988534365274  -  0,293912623038 

-  0,02681 12345019  -i-  0,109097496421  +  0,0016284675631 1 

-  0,0214631030242  +  0,0001456723382  U’® 

+  0,00167847752266  U" 

log  Ps  =  Y 


sc  <  0.0674:  Pf  =  10^ 

so  40:  Pr  =  0 

0.0674  <  sc  <  40: 

U  =  -  0,240657322658  +  1,36637719229  log  (sc) 

Y  =  +  3,40283217581  -  2,21030870597  U  -  0,218536586295 

+  0,895319589372  +  0,24989009775  -  0,569249436807 

-  0,1 1791682383  +  0,224131 16141 1  +  0,0245620259375 

-  0,04551 16002694  -  0,00190930738887 

+  0,00361471193389 

logPr  =  Y 


Igloos  en  zware  constructies: 

Ps  =  A’  *  (scp+CHsc)) 

Voor  A’,  B  en  C  gelden  de  waarden  als  in  tabel  D.l  gegeven. 


D.3 
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Tabel  D.l:  Curve  fit  coefficienten  voor  berekening  piekdruk  voor  igloos  en  zware  con- 
structies. 


voorzijde 

zijkant 

achterzijde 

A’ 

792,24 

227,84 

214,21 

B 

-2,0108 

-0,93655 

-1,1936 

C 

0,09419 

-0,12311 

-0,05778 

Impuls  (1.1-bronnen) 

Open  opslag  en  lichte  constructies: 
sc  <  0,0674:  Ig  =  10 

so  40:  Is  =  0 

0,0674  <  sc  <0,955: 

U  =  +  2,06761908721  +  3,0760329666  log  (sc) 

Y  =  +  2,52455620925  -  0,502992763686  U  +  0,171335645235 

+  0,0450176963051  -  0,01 18964626402 

0,955  <  sc  <  40: 

U  =  -  1,94708846747  +  2,40697745406  log  (sc) 

Y  =  +  1,67281645863  -  0,384519026965  U  -  0,0260816706301  U2 

+  0,00595798753822  +  0,014544526107  -  0,00663289334734 

-  0,00284189327204  +  0,013644816227 

log  Simp  =  Y 

Is  =  0.001  SimpQ^^3 

sc  <  0,0674:  If  =  100 

so  40:  Ir  =  0 

0,0674  <  sc  <  40: 

U  =  -  0,246208804814  +  1,33422049854  log  (sc) 

Y  =  +  2,70588058103  -  0,949516092853  U  +  0,112136118689  U2 

-  0,0250659183287 

log  Simp  =  Y 

Ir  =  0.001  Simp  Qi/3 


Igloos  en  zware  constructies: 

IJQ1I3  =A’  * 

Voor  A’,  B  en  C  gelden  de  waarden  als  in  tabel  D.2  gegeven. 

Tabel  D.2:  Curve  fit  coefficienten  voor  berekening  impuls  voor  igloos  en  zware  con¬ 
structies. 


voorzijde 

zijkant 

achterzijde 

A’ 

222,45 

122,90 

73,544 

B 

-0,87840 

-0,44293 

-0,36446 

C 

-0,03850 

-0,13156 

-0,12788 
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Schokgolfaankomsttijd  (1.1-bronnen) 

sc  <  0,674:  Ta  =  10-5 

so  40:  Ta  =  0 

0,0674  <  sc  <  40: 

U  =  -  0,202425716178  +  1,37784223635  log  (sc) 

Y  =  -  0,0591634288046  +  1,35706496258  U  +  0,052492798645  U2 

-  0,196563954086  -  0,0601770052288  +  0,0696360270891  U5 

+  0,0215297490092  -  0,0161658930785  -  0,00232531970294  U* 

+  0,00147752067524 U9 

log  Stim  =  Y 

Ta  =  0.001  StimQ‘/3 


Schokgolfsnelheid  (l.l-bronnen) 

sc  <  0,0674:  Vg  =  100 

so  40:  Vg  =  0 

0,0674  <  sc  <  40: 

U  =  -  0,202425716178  +  1,37784223635  log  (sc) 

Y  =  -  0,06621072854  -  0,698029762594  U  +  0,158916781906 

+  0,443812098136  -  0,1 13402023921  -  0,369887075049  U5 

+  0,129230567449  +  0,19857981 197  -  0,0867636217397 

-  0,0620391900135  +  0,0307482926566  +  0,0102657234407  1 

-  0,00546533250772  -  0,000693180974  3  +  0,0003847494916 

log  Svel  =  Y 

Vg  =  SvelQl^3 


Positieve  faseduur  (l.l-bronnen) 

sc  <  0,178:  Tp  =  10-4 

so  40:  Tp  =  0 

0,178  <  sc  <  1,01: 

U  =  +  1,92946154068  +  5,25099193925  log  (sc) 

Y  =  -  0,614227603559  +  0,130143717675  U  +  0,134872511954112 

+  0,0391574276906  -  0,00475933664702  -  0,00428144598008 

1,01  <  sc  <  2,78: 

U  =  -  2,12492525216  +  9,299628861 1  log  (sc) 

Y  =  +  0,315409245784  -  0,0297944268976  U  +  0,030632954288 

+  0,0183405574086  -  0,017396466621 1 -  0,00106321963633 

+  0,00562060030977  +  0,0001618217499  -  0,0006860188944 

2,78  <  sc  <  40: 

U  =  -  3,53626218091  +  3,46349745571  log  (sc) 

Y  =  +  0,686906642409  +  0,0933035304009  U  -  0,0005849420883  U2 

-  0,00226884995013  -  0,00295908591505  +  0,00148029868929 

log  Spos  =  Y 
Tp  =  0.001  SposQ*^3 
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HittestralingsefFect  (1.1  en  1.3-bronnen) 

Voor  1.3-explosieven  wordt  10%  van  de  totale  hoeveelheid  genomen,  terwijl  voor 
igloos  waarin  1.3  ligt  opgeslagen  de  hittestraling  alleen  wordt  berekend  aan  de 
voorzijde. 

Rc  =  2.8.QO-28 
Rg  =  0.45.Q0-^^ 
t  =  0.93.QO-2^ 

W=19.Q0-^/x2 

Indien  x  <  Rg  (object  in  de  vuurbal)  worden  de  waarden  voor  W  en  t  aangepast 
zodat  letaliteitsberekeningen  resulteren  in  kans  op  letaliteit  van  100%:  t  =  6; 

W  =  200. 


D.3  Formules  voor  (handmatige)  controle  van  de  resultaten  verkro’ 
gen  met  de  gefitte  formules 

Shock  velocity: 

V,  =  Co-V(l  +  6Ps  I  7Po) 

Particle  velocity: 

u  =  5.Ps  /  7Po.Vs  /  (1  +  6Ps/7Po) 

Dynamic  pressure: 

Pd  =  5/2Ps^l(7Po  +  Ps) 

Reflected  overpressure: 

Pr  =  2P,  [(7Po  +  4Ps)  /  (7Po  +  Ps)] 

Nu  zijn  alle  explosie-effecten  bekend  voor  deze  bron  (donor)  en  deze  afstand 
(object). 
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Bijlage  E  Formules  voor  de  berekening  van  de  ^Ballistic 
limit  velocities’ 


Voor  mensen  in  auto’s  wordt  aangenomen  dat  zij  worden  beschermd  door  een 
dunne  stalen  plaat  van  1,5  mm  dik.  Voor  de  penetratie  van  dunne  stalen  platen 
wordt  een  penetratiefunctie  gegeven  door  Healey  [11].  Als  deze  wordt  omgerekend 
naar  Sl-eenheden  en  voor  de  dikte  van  de  plaat  1,5  mm  wordt  ingevuld,  volgt  de 
formule  voor  de  ‘ballistic  limit  velocity’  voor  stalen  platen  welke  gebruikt  wordt  in 
het  programma: 

V50  =  41,11  (m/s) 

Dit  leidt  tot  ‘ballistic  limit  velocities’  van  76,5,  41,1  en  27,4  m/s  voor  penetratie 
door  staal  van  fragmenten  met  gewicht  respectievelijk  0,1, 1  en  4,5  kg. 

Voor  onbeschermde  personen  die  worden  getroffen  door  fragmenten  wordt  de 
‘ballistic  limit  velocity’  voor  huidperforatie  gebruikt.  Deze  formule  is  als  volgt 
[12]: 

V50  =  1247,1  Mm  +  22,03  (m/s) 
met 

m  =  massa  van  fragment  (kg) 

A  =  blootgestelde  oppervlakte  (m^) 

Er  wordt  aangenomen  dat  de  restsnelheid  na  huidperforatie  hoog  genoeg  is  om 
emstig  letsel  te  veroorzaken.  In  de  praktijk  is  het  niet  eenvoudig  om  een  waarde 
voor  A  te  definieren.  Daarom  wordt  een  veel  gebruikte  vormfactor  geintroduceerd: 
k  =  m/A3/2  [13].  Dit  leid  tot  de  volgende  formule: 

V50  =  1247,1  +  22,03  (m/s) 

Voor  de  meeste  effectieve  fragmenten  is  de  waarde  voor  k  gelijk  aan  4740  kg/m^^^. 
Voor  normaal  in  de  praktijk  voorkomende  fragmenten  is  deze  waarde  gelijk  aan 
2370  kg/ml''3.  Deze  laatste  waarde  wordt  gebruikt  in  het  programma. 
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Bijlage  F  Weergave  van  een  gedetailleerde  en  een  stan- 
daard  output  van  de  externe-risicoberekenin- 
gen 

Gebruikte  symbolen: 

P  =  invallende  druk  (kPa) 

I  =  invallende  impuls  (kPas) 

q  =  fragmentdichtheid  (l/m^) 

T  =  duur  van  de  positieve  fase  (s) 

Q  =  hittestraling  (kQ/m^) 

t  =  blootstellingstijd  tot  Q  (s) 

Ear  =  kans  op  oor-letsel  (%) 

Lung  =  kans  op  long-schade  (%) 

Frag  =  kans  op  letaliteit  door  fragmenten  (%) 

Frgm  =  kans  op  letaliteit  door  fragmenten  voor  fietsers/voetgangers  (%) 

Coll  =  kans  op  letaliteit  door  instorten  van  huizen  (%) 

Window  =  kans  op  letaliteit  door  ruitbreuk  (%) 

Body  =  kans  op  letaliteit  door  meesleureffect  als  gevolg  van  blast  (%) 

Bum  =  kans  op  letaliteit  door  hittestraling  (%) 

Bur2  =  kans  op  tweedegraads  brandwonden  (%) 

***  detailed  output 

Donor  IGLOOl  (Source  IGLOOl)  with  Contents: 

15000.  kg  1.1,  O.kg  1.2,  O.kg  1 .3  and  0.  kg  1 .4  (ihst=2) 

Overestimated  values  for  small  quantities  and  "not-open"  sites. 

Object  HSl  ,  Coordinates  (  75.,  600.),  Distance  464.,  Scaled  dist.  19. 

P=  6.5kPa,I=  415.Pa.s,q=  0.0  l/m2,  T=  144.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  54,  Body  0,  Bum  0,  Bur2  0 

Object  HS2  ,  Coordinates  (  150.,  690.),  Distance  553.,  Scaled  dist.  22. 

P=  5.2  kPa,  1=  350.  Pa.s,  q=  0.0  l/m2,  T=  15 1 .  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  12,  Body  0,  Bum  0,  Bur2  0 

Object  HS3  ,  Coordinates  (  230.,  600.),  Distance  476.,  Scaled  dist.  19. 

P=  6.3kPa,I=  405.Pa.s,q=  0.0  l/m2,  T=  145.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  47,  Body  0,  Bum  0,  Bur2  0 

Object  CM4  ,  Coordinates  (  550.,  200.),  Distance  437.,  Scaled  dist.  18. 

P=  7.1  kPa,  1=  440.  Pa.s,  q=  0.0  l/m2,  T=  141.  ms,  Q=  0.0  kW/ni2,  t=  0.0  s 
Prob:  Ear  l,Lung  0,  Frag  l,Frgm  OColl  0,  Window  0,  Body  0,  Bum  0,  Bur2  0 
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Object  HS6  ,  Coordinates  (  420.,  250.),  Distance  323.,  Scaled  dist.  13. 

P=  10.3  kPa,I=  589.  Pa.s,q=  0.1  l/ni2,  T=  129.  ms,  Q=  0.0  kW/ra2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  1,  Window  100,  Body  0,  Bum  0,  Bur2  0 


Object  HS8  ,  Coordinates  (  420.,  50.),  Distance  315.,  Scaled  dist.  13. 

P=  10.7  kPa,  1=  604.  Pa.s,  q=  0.1  l/m2,  T=  128.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  5,  Frgm  OColl  1,  Window  100,  Body  0,  Bum  0,  Bur2  0 


Object  HS9  ,  Coordinates  (  190.,  -245.),  Distance  389.,  Scaled  dist.  16. 

P=  8.2  kPa,  1=  492.  Pa.s,  q=  0.2  l/m2,  T=  136.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob;  Ear  0,  Lung  0,  Frag  0,  Frgm  OColl  0,  Window  91,  Body  0,  Bum  0,  Bur2  0 

Object  HS 10  ,  Coordinates  (  10., -245.),  Distance  397.,  Scaled  dist.  16. 

P=  7.9  kPa,  1=  482.  Pa.s,  q=  0.2  l/m2,  T=  137.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  89,  Body  0,  Bum  0,  Bur2  0 

Object  HSll  ,  Coordinates  ( -150.,  50.),  Distance  281.,  Scaled  dist.  11. 

P=  12.4  kPa,  1=  673.  Pa.s,  q=  0.2  l/m2,  T=  123.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  8,  Frgm  OColl  1,  Window  100,  Body  0,  Bum  0,  Bur2  0 

Object  RD2  ,  Coordinates  (  450.,  -100.),  Distance  409.,  Scaled  dist.  17. 

P=  7.7  kPa,  1=  469.  Pa.s,  q=  0.0  l/m2,  T=  138.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  l,Lung  0,  Frag  l,Frgm  2  Coll  0,  Window  32,  Body  0,  Bum  0,  Bur2  0 

Object  RD3  ,  Coordinates  (  395.,  -160.),  Distance  407.,  Scaled  dist.  16. 

P=  7.7  kPa,  1=  471.  Pa.s,  q=  0.0  l/m2,  T=  138.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  17,  Body  0,  Bum  0,  Bur2  0 


***  default  output 


Donor  IGLOO  1  and  source  IGLOO  1  give  lethality  or  injury  probabilities  for; 


Object  HSl  , 

type  HS, 

Object  HS2  , 

type  HS, 

Object  HS3  , 

type  HS, 

Object  CM4  , 

type  HU, 

Object  HS6  , 

type  HS, 

Object  HS8  , 

type  HS, 

Object  HS9  , 

type  HS, 

Object  HS  10, 

type  HS, 

Object  HSll  , 

type  HS, 

Object  RD2  , 

type  RD, 

Object  RD3  , 

type  RD, 

lethality;  0.03  and  *0.00, 

lethality;  0.01  and  0.00, 

lethality:  0.03  and  0.00, 

lethality:  0.01  and  0.00, 

lethality:  0.06  and  0.00, 

lethality:  0.10  and  0.00, 

lethality:  0.05  and  0.00, 

lethality;  0.05  and  0.00, 

lethality:  0.14  and  0.00, 

lethality;  0.03  and  0.02, 

lethality:  0.01  and  0.00, 


,  injury: 

0.06  and  *0.00 

injury: 

0.02  and  0.00 

injury; 

0.05  and  0.00 

injury: 

0.02  and  0.00 

injury: 

0.12  and  0.00 

injury: 

0.19  and  0.00 

injury: 

0.10  and  0.00 

injury: 

0.10  and  0.00 

injury: 

0.26  and  0.00 

injury: 

0.06  and  0.04 

injury: 

0.02  and  0.01 

*  second  figure  for  cyclists/pedestrians 
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Donor  IGLOO  1,  Type  13,  gives  7.0  Lethal  and  14.5  Serious  Injuries 
(Included  0.0  lethalities  on  roads) 

***  detailed  output 

Donor  EARTH2  (Source  EARTH2)  with  Contents: 

O.kg  1.1,  O.kg  1.2,  50000.  kg  1 .3  and  0.  kg  1 .4  (ihst=l) 

Overestimated  values  for  small  quantities  and  "not-open"  sites. 

***  default  output 

Donor  EARTH2  and  source  EARTH2  give  lethality  or  injury  probabilities  for: 


Donor  EARTH2,  Type  OB,  gives  0.0  Lethal  and  0.0  Serious  Injuries 
(Included  0.0  lethalities  on  roads) 

***  detailed  output 

Donor  IGL003  (Source  IGL003)  with  Contents: 

15000. kg  1.1,  O.kg  1.2,  oTkg  1 .3  and  0.  kg  1 .4 (ihst=2) 

Overestimated  values  for  small  quantities  and  "not-open"  sites. 

Object  HSl  ,  Coordinates  (  75.,  600.),  Distance  598.,  Scaled  dist.  24. 

P=  4.7  kPa,  1=  324.  Pa.s,  q=  0.0  l/m2,  T=  154.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  4,  Body  0,  Bum  0,  Bur2  0 

Object  HS3  ,  Coordinates  (  230.,  600.),  Distance  568.,  Scaled  dist.  23. 

P=  5.1  kPa,  1=  341.  Pa.s,  q=  0.0  l/m2,  T=  152.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  9,  Body  0,  Bum  0,  Bur2  0 

Object  CM4  ,  Coordinates  (  550.,  200.),  Distance  333.,  Scaled  dist.  13. 

P=  10.0  kPa,  1=  572.  Pa.s,  q=  0.0  l/m2,  T=  130.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  l,Lung  0,  Frag  2,  Frgm  OColl  0,  Window  0,  Body  0,  Bum  0,  Bur2  0 

Object  HS6  ,  Coordinates  (  420.,  250.),  Distance  269.,  Scaleddist.il. 

P=  13.2  kPa,I=  704.  Pa.s,q=  0.3  l/m2,  T=  121.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  11,  Frgm  OColl  1,  Window  100,  Body  0,  Bum  0,  Bur2  0 

Object  HS8  ,  Coordinates  (  420.,  50.),  Distance  158.,  Scaled  dist.  6. 

P=  28.4  kPa,I=  1163.  Pa.s,q=  0.7  l/m2,  T=  102.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  26,  Frgm  0  Coll  10,  WindowlOO,  Body  0,  Bum  0,  Bur2 
0 

Object  HS9  ,  Coordinates  (  190.,  -245.),  Distance  287.,  Scaled  dist.  12. 
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P=  12.1  kPa,  1=  661.  Pa.s,  q=  0.2  l/m2,  T=  124.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  7,  Frgm  OColl  1,  Window  100,  Body  0,  Bum  0,  Bur2  0 

Object  HS 10  ,  Coordinates  (  10., -245.),  Distance  375.,  Scaled  dist.  15. 

P=  8.5  kPa,  1=  510.  Pa.s,  q=  0.0  l/m2,  T=  135.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  95,  Body  0,  Bum  0,  Bur2  0 

Object  HSll  ,  Coordinates  ( -150.,  50.),  Distance  413.,  Scaled  dist.  17. 

P=  7.6kPa,I=  465.Pa.s,q=  0.2  l/m2,  T=  139.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  82,  Body  0,  Bum  0,  Bur2  0 

Object  RD2  ,  Coordinates  (  450., -100.),  Distance  230.,  Scaled  dist.  9. 

P=  16.4  kPa,  1=  819.  Pa.s,  q=  0.6  l/m2,  T=  1 15.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  l,Lung  0,  Frag  6,  Frgm  6  Coll  0,  Window  31,  Body  0,  Bum  0,  Bur2  0 

Object  RD3  ,  Coordinates  (  395., -160.),  Distance  234.,  Scaled  dist.  9. 

P=  16.0  kPa,I=  805.  Pa.s,q=  0.5  l/m2,  T=  116.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  6,  Frgm  6  Coll  0,  Window  20,  Body  0,  Bum  0,  Bur2  0 

***  default  output 


Donor  IGL003  and  source  IGL003  give  lethality  or  injury  probabilities  for: 


Object  HS3 

,  type  HS,  lethality: 

0.01  and  0.00,  injury: 

0.02  and  0.00 

Object  CM4 

,  type  HU,  lethality: 

0.02  and  0.00,  injury: 

0.03  and  0.00 

Object  HS6 

,  type  HS,  lethality: 

0.16  and  0.00,  injury: 

0.31  and  0.00 

Object  HS8 

,  type  HS,  lethality: 

0.37  and  0.00,  injury: 

0.65  and  0.00 

Object  HS9 

,  type  HS,  lethality: 

0.13  and  0.00,  injury: 

0.25  and  0.00 

Object  HSIO 

,  type  HS,  lethality: 

0.05  and  0.00,  injury: 

0.10  and  0.00 

Object  HSll 

,  type  HS,  lethality: 

0.05  and  0.00,  injury: 

0.09  and  0.00 

Object  RD2 

,  type  RD,  lethality: 

0.07  and  0.06,  injury: 

0.15  and  0.12 

Object  RD3 

,  type  RD,  lethality: 

0.07  and  0.06,  injury: 

0.13  and  0.12 

***  detailed  output 

Donor  IGL003  (Source  EARTH2)  with  Contents: 

0.  kg  1.1,  0.  kg  1.2,  50000.  kg  1.3  and  0.  kg  1.4  (ihst=l) 

Overestimated  values  for  small  quantities  and  "not-open"  sites. 


***  default  output 
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Donor  IGLOOS  and  source  EARTH2  give  lethality  or  injury  probabilities  for: 


***  detailed  output 

Donor  IGLOOS  (Source  OPEN4  )  with  Contents: 

0.  kg  1.1,  20000.  kg  1.2,  0.  kg  l.S  and  0.  kg  1.4  (ihst=l) 

Overestimated  values  for  small  quantities  and  "not-open"  sites. 

Object  RD2  ,  Coordinates  (  450., -100.),  Distance  200.,  Scaled  dist.  0. 

P=  0.0  kPa,  1=  0.  Pa.s,  q=  0.5  l/m2,  T=  0.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  6,  Frgm  6  Coll  0,  Window  0,  Body  0,  Bum  0,  Bur2  0 

Object  RDS  ,  Coordinates  (  S95., -160.),  Distance  186.,  Scaled  dist.  0. 

P=  0.0  kPa,  1=  0.  Pa.s,  q=  0.6  l/m2,  T=  0.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  7,  Frgm  7  Coll  0,  Window  0,  Body  0,  Bum  0,  Bur2  0 

***  default  output 

Donor  IGLOOS  and  source  OPEN4  give  lethality  or  injury  probabilities  for: 

Object  RD2  ,  type  RD,  lethality:  0.06  and  0.06,  injury:  0. IS  and  O.IS 

Object  RDS  ,  type  RD,  lethality:  0.07  and  0.07,  injury:  0.14  and  0.14 

Donor  IGLOOS,  Type  IS,  gives  16.7  Lethal  and  28. S  Serious  Injuries 
(Included  0.1  lethalities  on  roads) 

***  detailed  output 

Donor  OPEN4  (Source  OPEN4  )  with  Contents: 

O.kg  1.1,  20000.  kg  1.2,  0.  kg  l.S  and  0.  kg  1.4  (ihst=l) 

Overestimated  values  for  small  quantities  and  "not-open"  sites. 

Object  RD2  ,  Coordinates  (  450., -100.),  Distance  200.,  Scaled  dist.  0. 

P=  0.0  kPa,  1=  0.  Pa.s,  q=  0.5  l/m2,  T=  0.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 

Prob:  Ear  0,  Lung  0,  Frag  6,  Frgm  6  Coll  0,  Window  0,  Body  0,  Bum  0,  Bur2  0 

Object  RDS  ,  Coordinates  (  S95., -160.),  Distance  186.,  Scaled  dist.  0. 

P=  0.0  kPa,  1=  0.  Pa.s,  q=  0.6  l/m2,  T=  0.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 

Prob:  Ear  0,  Lung  0,  Frag  7,  Frgm  7  Coll  0,  Window  0,  Body  0,  Bum  0,  Bur2  0 

***  default  output 

Donor  OPEN4  and  source  OPEN4  give  lethality  or  injury  probabilities  for: 


Object  RD2  ,  type  RD,  lethality:  0.06  and  0.06,  injury:  O.IS  and  O.IS 


Object  RD3  ,  type  RD,  lethality:  0.07  and  0.07,  injury:  0.14  and  0.14 

Donor  OPEN4 ,  Type  OP,  gives  0.1  Lethal  and  0.1  Serious  Injuries 
(Included  0.0  lethalities  on  roads) 

***  detailed  output 

Donor  IGLOOS  (Source  IGLOOS)  with  Contents: 

ISOOO.  kg  1.1,  0.  kg  1.2,  0.  kg  1.3  and  0.  kg  1.4  (ihst=2) 

Overestimated  values  for  small  quantities  and  "not-open"  sites. 

Object  HS6  ,  Coordinates  (  420.,  2S0.),  Distance  433.,  Scaled  dist.  18. 

P=  7.1  kPa,  1=  444.  Pa.s,  q=  0.0  l/m2,  T=  141.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  72,  Body  0,  Bum  0,  Bur2  0 

Object  HS8  ,  Coordinates  (  420.,  SO.),  Distance  28S.,  Scaled  dist.  12. 

P=  12.2  kPa,  1=  66S.  Pa.s,  q=  0.2  l/m2,  T=  124.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  8,  Frgm  OColl  1,  Window  100,  Body  0,  Bum  0,  Bur2  0 

Object  HS9  ,  Coordinates  (  190., -24S.),  Distance  130.,  Scaled  dist.  S. 

P=  12.S  kPa,  1=  1393.  Pa.s,  q=  1.7  l/m2,  T=  9S.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  S3,  Frgm  0  Coll  1,  WindowlOO,  Body  0,  Bum  0,  Bur2  0 

Object  HS 10  ,  Coordinates  (  10., -24S.),  Distance  220.,  Scaled  dist.  9. 

P=  17.4  kPa,  1=  8S3.  Pa.s,  q=  0.6  l/m2,  T=  1 14.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  22,  Frgm  0  Coll  3,  WindowlOO,  Body  0,  Bum  0,  Bur2  0 

Object  HSU  ,  Coordinates  ( -ISO.,  SO.),  Distance  376.,  Scaled  dist.  IS. 

P=  8.S  kPa,  1=  S09.  Pa.s,  q=  0.3  l/m2,  T=  13S.  ms,  Q=  0.0kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  0,  Frgm  0  Coll  0,  Window  94,  Body  0,  Bum  0,  Bur2  0 

Object  RD2  ,  Coordinates  (  4S0., -100.),  Distance  263.,  Scaled  dist.  11. 

P=  13.6  kPa,I=  719.  Pa.s,q=  0.1  l/m2,  T=  121.  ms,  Q=  0.0kW/m2,t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  7,  Frgm  7  Coll  0,  Window  32,  Body  0,  Bum  0,  Bur2  0 

Object  RD3  ,  Coordinates  (  39S.,  ^160.),  Distance  212.,  Scaled  dist.  9. 

P=  18.3  kPa,  1=  883.  Pa.s,  q=  0.3  l/m2,  T=  1 12.  ms,  Q=  0.0  kW/m2,  t=  0.0  s 
Prob:  Ear  0,  Lung  0,  Frag  3,  Frgm  3  Coll  0,  Window  20,  Body  0,  Bum  0,  Bur2  0 

***  default  output 

Donor  IGLOOS  and  source  IGLOOS  give  lethality  or  injury  probabilities  for: 

Object  HS6  ,  type  HS,  lethality:  0.04  and  0.00,  injury:  0.08  and  0.00 

Object  HS8  ,  type  HS,  lethality:  0.13  and  0.00,  injury:  0.2S  and  0.00 
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Object  HS9  , 
Object  HS 10  , 
Object  HS 11  , 
Object  RD2  , 
Object  RD3  , 


type  HS,  lethality: 
type  HS,  lethality: 
type  HS,  lethality: 
type  RD,  lethality: 
type  RD,  lethality: 


0.56  and  0.00,  injury: 
0.29  and  0.00,  injury: 
0.05  and  0.00,  injury: 
0.09  and  0.07,  injury: 
0.04  and  0.03,  injury: 


1.00  and  0.00 
0.54  and  0.00 
0.11  and  0.00 
0.17  and  0.15 
0.08  and  0.06 


Donor  IGL005,  Type  13,  gives  9.2  Lethal  and  18.8  Serious  Injuries 
(Included  0.0  lethalities  on  roads) 
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Bijiage  G  Standaard  output  van  de  eindresultaten 
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He 
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Average  number  of  lethalities  for  this  complex  : 
Average  number  of  serious  injuries  : 

Collective  risk  : 


6.6 

12.3 

33.1  *  10-5 


* 

* 

* 
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Probability  of  lethality  (%) 


Code  of 
object 

Number  of 
objects 

Average 

Maximum 

Source  with 
maximum 

HSl 

5 

0.8 

3.0 

IGLOO  1 

HS3 

5 

0.8 

2.7 

IGLOOl 

CM4 

90 

0.7 

2.2 

IGL003 

HS6 

5 

5.2 

16.3 

IGL003 

HS8 

5 

12.0 

36.7 

IGL003 

HS9 

2 

14.9 

56.3 

IGL005 

HSIO 

1 

7.8 

28.7 

IGL005 

HSll 

1 

4.8 

13.7 

IGLOOl 

RD2 

100 

5.9 

12.7 

IGL003 

RD3 

10 

4.8 

12.8 

IGL003 
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